AN  INVESTIGATION  OF  CERTAIN 
THERMODYNAMIC  LOSSES  IN  MINIATURE 

CRYOCOOLERS 

FINAL  REPORT 


AUTHOR:  Cryogenics  Group,  Oxford  University 

RESEARCH:  Fast  Response  Temperature  Measurements  in  Cryocoolers 

REPORT  FOR:  USAF  AFRL 

CONTRACT:  FA8655-06-1-3071 

DATE:  02-May-2007 

ADDITIONAL  DOCUMENTS:  None. 

NUMBER  OF  PAGES :  7 1 

OUR  REF:  LOS-FR3 


REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  Information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply 
with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  ('DD-M/W-yyvr;  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  Toj 

08-05-2007  Final  Report  21  July  2006  -  21  -Mar-07 


4.  TITLE  AND  SUBTITLE 


5a.  CONTRACT  NUMBER 

FA8655-06-1-3071 


5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 


5d.  PROJECT  NUMBER 


Dr.  Charles  R.  Stone 


5d.  TASK  NUMBER 


5e.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Oxford 
Parks  Road 
Oxford  0X1  3PJ 
United  Kingdom 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


EOARD 

Unit  4515  BOX  14 
APO  AE  09421 


11.  SPONSOR/MONITOR’S  REPORT  NUMBER(S) 

Grant  06-3071 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

Measurements  were  made  of  gas  temperatures  and  power  losses  for  a  compression  volume  that  included  a  regenerative  heat  exchanger. 

A  fine  wire  resistance  thermometer  with  a  rapid  response  was  used  to  measure  the  gas  temperature  for  flows  between  the  piston/cylinder 
volume  and  the  regenerator.  The  temperature  profiles  exhibited  features  that  varied  with  different  gases  and  operating  conditions.  These 
results  were  interpreted  in  terms  of  a  stratified  tidal  gas  volume  that  is  subject  to  thermal  diffusion  and  flow  turbulence  effects.  A  few  nitrogen 
results  for  high  pressure  and  low  frequency  were  very  unstable  with  significant  high  frequency  components.  A  tentative  interpretation  was  that 
this  behaviour  was  associated  with  the  transition  from  laminar  to  turbulent  flow  conditions. 

At  the  other  end  of  the  regenerator  in  a  separate  fixed  volume,  a  twin  thermocouple  probe  was  used  to  investigate  the  performance  of  a 
compensation  technique  designed  to  improve  the  time  response  of  the  bare  thermocouples.  The  technique  involves  arranging  two 
thermocouples  with  different  time  constants  such  that  they  experience  the  same  gas  temperature  variation.  The  gas  temperature  is 
reconstructed  by  using  various  approaches  for  estimating  the  response  times.  It  is  shown  that  this  technique  can  be  considered,  provided  that 
the  noise  levels  are  low  enough.  For  the  measurements  described  here  the  noise  could  only  be  sufficiently  reduced  by  ensemble  averaging. 
This  would  be  a  limitation  in  circumstances  where  the  waveform  is  not  very  repeatable. 

The  power  loss  measurements  showed  that  a  loss  existed  that  could  be  correlated  by  a  combination  of  Reynolds  number  and  the 
experimental  parameters  pressure,  frequency  and  stroke.  Attempts  to  formulate  this  in  a  non-dimensional  form  were  not  successful 
suggesting  that  key  parameters  have  been  missed.  One  aspect  in  particular  that  needs  closer  attention  is  the  effect  of  phase  angle  between 
mass  flow  and  pressure  variation. 

An  interesting  feature  of  this  work  is  the  general  agreement  between  the  resistance  wire  thermometry  results  and  the  power  loss 
measurements.  Both  sets  of  results  exhibit  changes  in  behaviour  consistent  with  the  idea  of  a  laminar  to  turbulent  transition.  Furthermore  a 
region  of  transition  defined  by  a  Reynolds  number  suggested  by  the  power  measurements  is  reasonably  supported  by  the  thermometer 
results.. 


15.  SUBJECT  TERMS 

Cryocoolers,  Thermodynamics,  Compressible  Flows,  EOARD 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

UNCLAS 

b.  ABSTRACT 

UNCLAS 

c.  THIS  PAGE 

UNCLAS 

17.  LIMITATION  OF 
ABSTRACT 

UL 


18,  NUMBER 
OF  PAGES 

72 


19a.  NAME  OF  RESPONSIBLE  PERSON 

BARRETT  A.  FLAKE 

19b.  TELEPHONE  NUMBER  (include  area  code) 

+44  (0)1895  616144 


Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39-18 


1  CONTENTS 


1  CONTENTS . 2 

2  EXECUTIVE  SUMMARY . 3 

3  BACKGROUND . 4 

4  INTRODUCTION . 4 

5  EXPERIMENT AE  SETUP . 7 

6  PINE  WIRE  RESISTANCE  THERMOMETER  MEASUREMENTS  (HOT  WIRE 

PROBE) . 10 

7  MEASUREMENTS  MADE  USING  COMPENSATED  THERMOCOUPEES . 19 

8  DESCRIPTION  OP  POWER  LOSS  MEASUREMENTS . 22 

9  OVERALL  SUMMARY  AND  CONCLUSIONS . 37 

10  FUTURE  WORK . 38 

11  REFERENCES . 38 

APPENDIX  I:  PRESSURE  DROP  IN  REGENERATOR . 39 

APPENDIX  2:  REPORT  ON  COMPENSATED  THERMOCOUPLE  RESULTS . 42 


2  EXECUTIVE  SUMMARY 

Measurements  were  made  of  gas  temperatures  and  power  losses  for  a  compression  volume 
that  included  a  regenerative  heat  exchanger. 

A  fine  wire  resistance  thermometer  with  a  rapid  response  was  used  to  measure  the  gas 
temperature  for  flows  between  the  piston/cy finder  volume  and  the  regenerator.  The 
temperature  profiles  exhibited  features  that  varied  with  different  gases  and  operating 
conditions.  These  results  were  interpreted  in  terms  of  a  stratified  tidal  gas  volume  that  is 
subject  to  thermal  diffusion  and  flow  turbulence  effects.  A  few  nitrogen  results  for  high 
pressure  and  low  frequency  were  very  unstable  with  significant  high  frequency  components. 

A  tentative  interpretation  was  that  this  behaviour  was  associated  with  the  transition  from 
laminar  to  turbulent  flow  conditions. 

At  the  other  end  of  the  regenerator  in  a  separate  fixed  volume,  a  twin  thermocouple  probe  was 
used  to  investigate  the  performance  of  a  compensation  technique  designed  to  improve  the 
time  response  of  the  bare  thermocouples.  The  technique  involves  arranging  two 
thermocouples  with  different  time  constants  such  that  they  experience  the  same  gas 
temperature  variation.  The  gas  temperature  is  reconstructed  by  using  various  approaches  for 
estimating  the  response  times.  It  is  shown  that  this  technique  can  be  considered,  provided  that 
the  noise  levels  are  low  enough.  For  the  measurements  described  here  the  noise  could  only  be 
sufficiently  reduced  by  ensemble  averaging.  This  would  be  a  limitation  in  circumstances 
where  the  waveform  is  not  very  repeatable. 

The  power  loss  measurements  showed  that  a  loss  existed  that  could  be  correlated  by  a 
combination  of  Reynolds  number  and  the  experimental  parameters  pressure,  frequency  and 
stroke.  Attempts  to  formulate  this  in  a  non-dimensional  form  were  not  successful  suggesting 
that  key  parameters  have  been  missed.  One  aspect  in  particular  that  needs  closer  attention  is 
the  effect  of  phase  angle  between  mass  flow  and  pressure  variation. 

An  interesting  feature  of  this  work  is  the  general  agreement  between  the  resistance  wire 
thermometry  results  and  the  power  loss  measurements.  Both  sets  of  results  exhibit  changes  in 
behaviour  consistent  with  the  idea  of  a  laminar  to  turbulent  transition.  Furthermore  a  region  of 
transition  defined  by  a  Reynolds  number  suggested  by  the  power  measurements  is  reasonably 
supported  by  the  thermometer  results. 


3  BACKGROUND 


This  report  describes  work  that  is  a  continuation  of  previous  studies  reported  in  Los-FRl  ^  and 
Los-FR2^.  These  investigations  are  intended  to  give  a  better  understanding  of  thermodynamic 
losses  that  occur  in  miniature  Stirling  cycle  coolers  but  which  are  not  adequately  explained  by 
our  existing  models.  It  is  hoped  that  the  results  of  this  study  will  be  useful  in  two  ways.  Firstly 
the  results  will  help  to  develop  improved  CFD  models  of  processes  that  occur  in  these 
coolers;  secondly  that  the  results  might  produce  ideas  for  reducing  such  losses. 

Los-FRl  looked  closely  at  methods  for  determining  power  input  to  the  gas.  The  conclusion 
drawn  was  that  the  only  reliable  way  was  to  determine  the  integral  of  P.dV  for  the 
compression  space.  This  approach  is  accurate  providing  that  any  phase  changes  introduced  in 
the  instrumentation  electronics  are  corrected.  It  was  found  that  there  was  too  much 
uncertainty  in  the  motor  losses  to  allow  the  power  input  to  the  gas  to  be  determined  from  the 
total  input  power. 

Los-FRZ  used  the  test  rig  configured  as  a  simple  gas  spring  to  produce  gas  springs  losses  that 
could  be  compared  with  established  results,  in  particular  the  work  of  Kornhauser^’"^.  It  was 
found  that  the  agreement  was  good  for  a  wide  range  of  gases  and  Peclet  numbers,  and  that  the 
results  supported  Kornhauser’s  modification  to  Lee’s^  original  equation.  For  these 
measurements  there  was  a  significant  gas  leakage  past  the  piston  and  the  results  also 
demonstrated  that  this  seal  leakage  did  not  appear  to  affect  the  nature  of  the  gas  spring  losses. 
LOS-FR2  also  briefly  described  an  unsuccessful  attempt  to  apply  a  compensated  twin 
thermocouple  sensor  to  the  measurement  of  varying  gas  temperature. 

The  final  part  of  LOS-FR2  looked  at  whether  the  measurement  of  compression  losses  could 
be  used  to  deduce  heat  transfer  coefficients  in  a  regenerator  for  oscillating  flow.  A  differential 
method  was  suggested  but  not  implemented. 

In  the  work  described  here,  the  areas  of  temperature  measurement  and  gas  spring  power  loss 
are  again  pursued,  this  time  with  a  gas  spring  geometry  that  more  closely  resembles  the 
geometry  of  a  real  Stirling  machine. 


4  INTRODUCTION 

The  understanding  of  the  compression  losses  in  Stirling  cycle  machines  can  be  approached  in 
two  different  ways: 

•  A  Quantitive  approach:  Direct  measurement  of  losses  -  Losses  can  be  analysed  to 
derive  values  for  components  that  are  not  included  in  existing  models.  The  behaviour 
of  these  residual  losses  can  be  investigated  for  varying  compressor  parameters 

•  A  Qualitative  approach:  Measurement  of  cycle  temperature  variation  at  particular 
locations.  Temperature  measurements  can  indicate  magnitudes  of  heat  transfer  and 
mixing  processes. 


These  are  not  in  principle  entirely  separate  -  accurate  temperature  measurements  for  enough 
points  would  allow  the  energy  flows  and  hence  losses  to  be  determined.  However,  in  reality  it 
is  not  practical  to  do  this  -  temperature  sensors  tend  to  be  intrusive  and  only  a  limited  number 
can  be  installed. 

The  emphasis  that  will  be  pursued  here  is  a  closer  consideration  of  the  problems  involved 
with  measuring  varying  gas  temperatures  in  Stirling  cycle  machines.  Two  approaches  were 
investigated:  a  fine  wire  resistance  thermometer  (RT)  and  a  more  sophisticated  technique  that 
will  be  referred  to  as  a  compensated  thermocouple  thermometer. 

To  provide  a  suitable  test  bed  for  the  thermometry,  the  test  rig  was  configured  as  a  more 
complex  gas  spring  that  included  a  regenerative  heat  exchanger.  Additional  components  were 
made  to  allow  the  different  thermometers  to  be  installed.  The  modified  test  rig  is  shown  in 
Figure  1. 


Figure  1,  Compression  Test  Rig  with  New  Configuration 


Although  the  main  focus  of  this  work  was  the  investigation  of  thermometry  it  was  decided 
that  it  would  also  be  worth  measuring  the  gas  spring  losses  for  this  more  complicated 
configuration.  The  previous  work  on  the  simple  gas  spring  demonstrated  that  the  phase 
correction  approach  can  give  good  results  and  this  was  an  opportunity  to  obtain  some 
additional  data. 

To  summarize  it  is  hoped  that  although  information  acquired  with  these  measurements  will 
not  be  enough  to  guarantee  a  clear  description,  they  will  nonetheless  be  useful  in  gaining  a 
better  understanding: 

•  Temperature  measurements  are  point  measurements  so  they  cannot  be  used  to 

determine  energy  flows  without  assuming  particular  temperature  profiles.  For  these 


oscillating  flow  conditions  and  the  small  mean  differenees  that  would  have  to  be 
determined,  assuming  a  temperature  distribution  is  not  justifiable  -  uncertainties 
would  be  too  large  for  the  results  to  be  useful 

•  The  cold  wire  RT  has  a  fast  response  and  will  give  temperature  variation  for  a  small 
volume  near  the  regenerator.  This  will  be  useful  for  investigating  the  nature  of  the 
flows  between  the  compression  space  and  the  regenerator. 

•  The  twin  thermoeouple  is  expeeted  to  have  insuffieient  response  to  follow  temperature 
variation.  However  it  is  a  useful  test  vehiele  for  the  compensated  thermocouple 
approach.  It  is  not  expected  that  it  will  give  as  much  insight  into  oscillating  flows  as 
the  cold  wire  RT  beeause  the  flow  here  will  be  small.  Nonetheless  it  will  give 
experience  that  will  help  determine  whether  this  teehnique  can  be  applied  to  these 
types  of  measurements 

•  Power  loss  measurements  will  help  understand  whether  losses  ean  be  eharaeterised  in 
a  useful  way.  They  will  not  necessarily  be  directly  applicable  to  Stirling  eycle 
machines  but  they  will  be  eonsiderably  more  realistic  than  those  calculated  on  the 
basis  of  a  simple  gas  spring  model. 


Old  Configuration 
Simple  Gas  Spring 


New  Configuration 
hcludes  Regenerator 


Figure  2,  Comparison  between  Old  and  New  Configurations 

A  cross  section  of  the  test  rig  used  in  this  work  is  shown  in  section,  together  with  the  original 
simple  gas  spring  in  Figure  2.  It  will  be  seen  that  for  the  new  configuration  the  compression 
spaee  above  the  piston  has  been  redueed  in  volume  and  is  now  attaehed  to  two  additional 


cylindrical  volumes.  The  one  elosest  to  the  eompression  spaee  is  filled  with  metal  mesh  so 
that  it  functions  as  a  regenerative  heat  exchanger  (i.e.  a  regenerator),  the  seeond  added  volume 
is  free  of  any  additional  heat  transfer  surface. 

The  fine  wire  resistanee  thermometer  was  positioned  slightly  off  eentre  in  the  top  of  the 
eompression  spaee  (Volume  A)  elose  to  the  regenerator  (Volume  B).  The  temperature  that 
will  be  measured  will  be  a  function  of  both  the  heat  transfer  eonditions  in  the  eompression 
volume  and  the  temperature  of  the  gas  flowing  into  the  volume  from  the  regenerator. 

The  eompensated  thermoeouple  thermometer  was  positioned  approximately  in  the  eentre  of 
the  top  volume  (Volume  C).  The  gas  flow  into  the  volume  is  very  mueh  smaller  than  for  the 
volume  A  and  it  is  expeeted  that  the  temperature  variation  here  will  be  less  affected  by  the 
flow.  Neither  of  these  thermometers  was  purpose  built  for  this  applieation.  The  positioning 
was  therefore  mainly  dietated  by  their  existing  installation  requirements. 

The  power  loss  measurements  were  made  by  reeording  the  instantaneous  values  of  piston 
position  and  eompression  spaee  pressure.  After  eorreeting  for  phase  ehanges  in  the  eleetronies 
the  power  loss  is  ealeulated  from: 


Power  =  /.(^  P.dV 

The  volume  ehanges  are  ealeulated  from  the  piston  position  and  piston  area 

Temperature  and  power  loss  measurements  were  made  for  helium  and  nitrogen  for  varying 
pressures,  frequeneies  and  strokes.  A  Fylde  8  channel  data  aequisition  system  was  used  to 
eollect  the  data.  This  will  be  deseribed  in  more  detail  in  the  next  seetion. 


5  EXPERIMENTAL  SETUP 

A  detailed  diagram  of  the  test  rig  volumes  with  dimensions  is  show  in  Figure  3.  The  prineipal 
additions  are  the  regenerator  volume  and  the  two  thermometer  installations. 

Regenerator 

The  regenerator  volume  was  paeked  with  eireular  mesh  discs  punched  out  from  mesh  sheet. 
For  the  measurements  deseribed  only  one  type  of  mesh  was  used  (see  Table  1)  but  the  rig 
eould  be  repaeked  with  other  types  of  mesh  if  required. 

Table  1. 


Regenerator 

Material 

Mesh  Weave 

Mesh  Number 
(wires  per 
inch) 

Wire  Diameter 
(mm) 

Porosity* 

(void  vol/total  vol) 

Stainless  Steel 

Plain  Weave 

#250 

0.040 

0.67 

*The  porosity  was  determined  by  weighing  the  regenerator  assembly  before  and  after  packing 
The  mass  of  regenerator  material  and  its  density  gives  the  volume  of  regenerator  mesh. 


Handling  loose  mesh  dises  ean  be  problematie  beeause  the  material  as  woven  is  very  springy 
To  help  with  paeking  the  regenerator  the  mesh  was  first  annealed  in  a  vaeuum  furnaee.  The 
annealed  mesh  was  then  paeked  in  an  assembly  jig  prior  to  being  installed  in  the  regenerator 


housing.  To  prevent  the  regenerator  from  bowing  out  at  its  ends  two  supports  were  included 
as  shown  in  Figure  4.  The  supports  are  designed  to  be  as  open  as  possible  so  as  to  avoid 
affecting  the  gas  flows  into  and  out  of  the  regenerator.  As  it  turned  it  was  found  that  the 
annealing  process  caused  a  slight  bonding  of  the  meshes  together  so  that  they  handled  as 
slightly  sintered  blocks.  The  supports  were  therefore  probably  not  really  necessary. 


Figure  3,  Details  of  New  Configuration  for  Test  Rig 


Area  1 
Perim 

Total 

Perim 

Area2 

Perim 


=  24.37 
eterl  =  1  9.00 

Area=  1  78.41 
eter=  1  34.01 


Figure  4,  Regenerator  Mesh  Support 


Thermometer  installations 

The  thermometers  required  speeial  adaptors  to  be  made  to  allow  them  to  be  installed  with  the 
minimum  of  extra  dead  volume  and  intrusion  into  the  gas  flow.  These  are  shown  in  figure  7. 

Data  Aquisition  System 

The  data  aequisition  equipment  used  was  the  Fylde  8  ehannel  unit  used  previously.  The 
output  of  the  Fylde  unit  is  direetly  imported  in  to  a  Excel  spread  sheet  for  subsequent 
analysis.  Because  it  was  not  intended  to  use  the  electrical  input  power  for  loss  calculations, 
the  instantaneous  values  for  motor  current  and  voltage  measurements  were  no  longer 
recorded.  Additional  inputs  to  the  A/D  were  the  thermometer  outputs  and  an  extra  pressure 
channel  to  record  the  pressure  variation  at  the  end  of  the  regenerator.  The  channel  allocation 
with  sensor  details  and  calibrations  are  show  in  Table  2. 

Table  2:  Channel  allocations  for  inputs  into  Fylde  Data  Acquisition  system 


Channel 

Measurement  Parameter 

Sensor 

Gain 

Offset 

1 

Pressure  in  Volume  A 

Endevco 

2.408  bar/V 

6.0206  bar 

2 

Pressure  in  Volume  C 

Druck  200 

1.9677  bar/V 

15.444  bar 

3 

N/C 

4 

Pressure  in  compressor  body 

Druck  820 

6.0046  bar/V 

8.3174  bar 

5 

Piston  Position 

Sch’  EVDT 

2.0884  mm/V 

-5.5968  mm 

6 

Fine  Wire  Thermometer 

TSI  probe 

10.296  K/V 

38.379  degC 

7 

TC  25  micron 

Custom  probe 

4.00  K/V 

24.00  deg  C 

8 

TC  50  micron 

Custom  probe 

4.00  K/V 

24.00  deg  C 

Note:  The  Fylde  data  Acquisition  unit  digitises  two  channels  at  a  time  with  5  microsecond 
intervals  hence  channels  1  and  5  are  at  t=0,  channels  2  and  6  are  at  t=5  microseconds  etc 


In  addition  to  the  cyclic  variations  recorded  by  the  Fylde  unit,  information  was  also  recorded 
manually: 

Gas  used 

Frequency  (from  power  analyser) 

Total  Input  Power 
Coil  Current 
Nominal  Fill  Pressure 
Regenerator  Body  Temperature 
Cold  Junction  Temperature  (for  TCs) 

Amp7  zeros  (V)  TC50 
Amp8  zeros  (V)  TC50 
Resistance  Thermometer  current 

Experimental  Procedure 

The  procedure  followed  was  to  fdl  the  test  rig  with  the  chosen  gas  to  a  prescribed  pressure. 
The  frequency  was  set  and  the  stroke  was  increased  to  the  required  value  by  observing  the 
LVDT  output  on  a  scope.  The  mean  position  of  the  piston  was  maintained  approximately  in 
its  mid  position  by  varying  the  dc  current  through  the  motor  coil.  When  the  compressor 
operation  reached  a  steady  state  then  the  Fylde  data  acquisition  system  was  used  to  record  the 
variation  of  the  parameters  defined  in  the  channel  allocation  above.  The  manual 
measurements  were  also  recorded. 

Two  different  setups  were  used.  In  one,  the  sampling  period  was  fixed  at  1  second  with  a 
sample  rate  of  5kHz.  These  measurements  were  used  for  the  resistance  thermometry  and 
power  measurements.  For  the  compensated  thermocouple  measurements  it  was  desirable  to 
enable  greater  ensemble  averaging.  The  sample  rate  was  reduced  to  2kHz  and  the  sample 
period  was  increased  to  15  seconds. 

The  measurements  recorded  in  the  Excel  spreadsheet  were  the  voltage  outputs  taken  directly 
from  the  various  instruments.  Phase  compensation  and  calibration  factors  are  applied  to  these 
values  to  arrive  at  converted  parameter  values  that  can  be  used  for  analysis.  This  process  of 
phase  compensation  and  calibration  is  executed  in  additional  Excel  spreadsheets. 


6  FINE  WIRE  RESISTANCE  THERMOMETER  MEASUREMENTS 
(HOT  WIRE  PROBE) 

Description  of  Probe 

The  resistance  thermometer  measurements  were  made  using  a  standard  hot  wire  anemometer 
probe  manufactured  by  TSI.  The  complete  probe  assembly  consists  of  two  parts,  a  sensor  that 
incorporates  the  fine  wire  sensing  element  and  a  support  into  which  the  sensor  component 
fits.  The  manufacturer’s  model  numbers  and  details  are  as  follows: 

Hot  Wire  Sensor:  TSI  Model  1210-T1.5  General  Purpose  Probe  with  wire  sensor.  (The 

wire  is  3.8  microns  in  diameter  and  is  tungsten  with  a  protective 
platinum  coating) 

Probe  Support:  TSI  Model  1150  Single  Sensor,  Standard  Probe  Support 

(The  probe  supports  come  in  three  lengths,  the  model  used  was 
standard  length  =  6  inch  version  designated  Model  1150-6) 


Hot  wire  sensor  and  probe  support  are  shown  in  Figures  5  and  6  respeetively 


© 


Figure  5,  TSI  Model  1210-T1.5  General  Purpose  Probe  With  Wire  Sensor. 


Figure  6,  TSI  Model  1150  Single  Sensor,  Standard  Probe  Support 

For  use  as  an  anemometer  this  type  of  probe  is  usually  operated  with  a  varying  eurrent  that 
maintains  a  eonstant  temperature.  However  to  operate  it  as  a  thermometer  the  sensor  is 
energized  with  a  eonstant  eurrent.  For  the  results  deseribed  here  the  eurrent  was  1mA.  The 
ehoiee  of  eurrent  is  a  eompromise  between  sensitivity  and  errors  produeed  by  self  heating.  For 
the  nitrogen  at  1  bar,  the  self  heating  was  found  to  be  aeeeptably  small  by  observing  that  there 
was  no  obvious  ehange  in  temperature  variation  when  the  eurrent  was  redueed  to  0.5mA. 

The  output  from  the  probe  eonditioning  eleetronies  is  the  voltage  aeross  the  wire  probe. 

Henee  if  the  temperature/resistanee  eharaeteristies  of  the  probe  are  known  then  the 
temperature  ean  be  ealeulated.  The  nominal  resistanee  is  6  ohms  and  the  sensitivity  with  1 
mA  is  ~  10.3  deg  C  per  volt.  This  output  was  reeorded  using  the  Fylde  Data  Aequisition  unit 
The  response  time  for  this  thermometer  is  low  beeause  of  the  very  small  sensor  diameter.  It 
was  intended  to  estimate  this  by  observing  the  temperature  deeay  when  the  sensor  is 
temporarily  heated  by  a  short  pulse.  Unfortunately  the  faeility  for  doing  this  in  the  TSI 
instrument  did  not  appear  to  work.  Therefore  the  time  eonstant  is  estimated  using  a  heat 
transfer  eorrelation  for  a  eylinder.  With  this  small  diameter,  the  Reynolds  number  is  mueh 
smaller  than  1 ,  so  the  flow  is  laminar  and  unseparated.  The  Nusselt  number  is  evaluated  using 
Churehill  and  Bernstein’s  equation: 


0.62Re;^"Pr''' 

l  +  (0.4/Pr)"''J'^ 

Sinee  the  temperature  fluetuations  are  small,  the  temperature  sensitivity  of  fluid  properties  are 
not  taken  into  aeeount.  All  fluid  properties  are  evaluated  at  25  deg  C.  From  the  Nusselt 
number  definition,  the  heat  transfer  eoeffieient  Qi)  ean  be  derived,  whieh  is  then  used  to 
ealeulate  the  time  eonstant  as  shown  below. 
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where  the  density  (/y,)  and  speeifie  heat  (Cw)  are  values  for  tungsten. 

Figure  7  shows  that  the  time  eonstant  reduees  with  inereasing  pressure.  Due  to  the  mueh 
higher  thermal  eonduetivity  of  helium,  the  eold  wire  has  a  slower  response  in  a  nitrogen 
environment.  Therefore,  the  biggest  time  eonstant  oeeurs  at  1  bar  pressure  in  nitrogen  gas  and 
zero  flow  speed.  The  time  eonstant  is  still  1 .2  ms,  whieh  is  small  enough  for  temperature 
measurements  to  be  eonsidered  instantaneous  at  the  operating  frequeneies  used  in  this  study 
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Figure  7,  Estimated  time  constant  of  cold  wire  at  0,5  m/s  flow  velocity. 

The  eomplete  probe  assembly  was  installed  into  the  eompression  test  rig  using  a  redesigned 
eompression  spaee  plate  and  a  housing  that  provided  a  seal  on  the  OD  of  the  4.6  mm  probe 
support.  The  housing  was  designed  to  allow  the  probe  to  be  positioned  in  different  radial 
loeations  within  the  eompression  volume  but  for  the  results  given  here  the  probe  was 
positioned  radially  at  the  eentre  and  axially  about  4.5  mm  away  from  the  regenerator.  The 
installation  of  the  probe  and  probe  housing  are  shown  in  Figure  8. 
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Figure  8,  Installation  Details  for  Pressure  Transducers  and  Thermometers 


Gas  Temperature  Measurements 

The  gas  temperature  variation  was  measured  over  a  wide  range  of  conditions  for  the  two 
gases,  nitrogen  and  helium.  The  ranges  covered  were; 

Pressure;  1  to  26  bar 

Frequency;  8.6  to  53  Hz 

Stroke;  4  to  14  mm 

It  was  not  possible  to  obtain  complete  coverage  of  the  strokes  for  all  pressure  and  frequency 
combinations  because  the  motor  was  largely  running  off  resonance  and  the  coil  current  had  to 
be  limited  to  prevent  the  coil  overheating. 

Overall  the  measurements  revealed  a  surprisingly  complex  variation  in  temperature  profiles. 
Part  of  the  variation  that  is  seen  with  stroke  is  likely  to  be  due  to  the  varying  time  that  the 
sensor  spends  in  the  tidal  gas  flow  that  is  flowing  in  and  out  of  the  regenerator.  As  the  sensor 
is  approximately  4.5  mm  away  from  the  regenerator,  the  sensor  would  be  expected  to  respond 
to  the  tidal  flow  only  for  small  periods  if  the  stroke  is  small  (e.g.  6  mm).  When  the  stroke 
rises  to  14  mm  then  the  tidal  flow  will  dominate  the  temperature  profile. 


It  has  not  been  possible  to  provide  a  comprehensive  explanation  for  all  that  has  been 
observed.  There  is  some  overlap  between  the  nitrogen  and  helium  temperature  profiles  but 
there  are  also  regimes  where  the  nitrogen  and  helium  data  were  found  to  be  quite  distinct.  The 
most  general  observation  that  was  common  to  both  sets  of  data,  was  that  whilst  cycle  to  cycle 
repeatability  for  the  temperature  could  be  very  variable  this  was  not  reflected  in  the  pressure 
or  piston  displacement  values  which  were  very  repeatable  and  generally  sinusoidal. 

Helium  Results 

The  range  of  temperature  profiles  obtain  with  helium  is  fairly  well  illustrated  in  Figures  9  and 
10  which  show  the  variation  for  differing  strokes.  The  results  in  Figure  9  represent  the  data 
with  the  lowest  mass  flows  i.e.  low  frequency  and  pressure.  It  is  seen  that  whilst  the 
waveform  becomes  more  distorted  as  the  stroke  increases  up  to  the  maximum  of  14mm,  there 
are  no  additional  features  in  the  waveform.  It  is  possible  however  to  detect  discontinuities  in 
the  gradient  for  the  14  mm  stroke. 


8.5  Hz,  1  bar  pressure 


Figure  9,  Fine  Wire  RT  measurements:  Helium,  Ibar  8,5Hz  for  different  strokes  (mm) 


In  Figure  10  the  results  are  shown  for  a  pressure  of  5  bar,  still  with  a  low  frequency  of  8.5  Hz. 
The  waveform  has  acquired  the  feature  of  an  extra  minimum  in  the  low  temperature  end.  At 
higher  frequencies  the  temperature  profile  changes  but  it  still  retains  the  same  general  features 
as  is  shown  in  Figure  1 1 . 

At  the  highest  pressures  and  strokes  the  waveform  starts  to  show  a  distinct  lack  of 
repeatability  with  apparent  high  frequency  components  as  is  shown  in  Figure  12.  It  is  seen 
that  the  second  peak  is  quite  abrupt  and  that  there  are  signs  of  instability. 


Number  of  cycles 


Figure  10,  Fine  Wire  RT  measurements:  Helium,  5bar,  8.5Hz  for  different  strokes 

(mm) 


Helium 


Figure  11,  Fine  Wire  RT  measurements:  Helium,  3bar,  29  Hz,  stroke  =  14mm 
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Figure  12,  Fine  Wire  RT  measurements:  Helium,  21bar,  44  Hz,  stroke  =  14mm 


Essentially  the  behaviour  can  be  summarised  as  follows; 

•  At  the  lowest  pressures  the  temperature  variation  is  basically  sinusoidal  with  distortion 
increasing  with  stroke  and  some  sign  of  a  discontinuity  in  the  gradient 

•  As  the  pressure  is  increased,  the  waveform  distortion  is  also  increased  until  a  second 
minimum  appears  in  the  lower  temperature  region.  Discontinuities  in  the  gradient  also 
become  more  obvious. 

•  The  results  do  not  allow  the  transition  to  be  pin-pointed  because  the  next  pressure 
investigated  was  3  bar  and  the  double  minima  appear  to  be  established  for  14  mm. 

•  As  the  frequency  and  pressure  is  increased  the  second  minimum  remains  but  its  shape 
and  position  change 

•  At  the  highest  pressures,  frequencies  and  strokes  the  second  minima  becomes 
increasingly  abrupt.  Also  the  waveform  show  signs  of  instability  with  high  frequency 
components  becoming  apparent. 


Nitrogen  Results 

The  low  pressure  nitrogen  results  are  found  to  be  fairly  similar  to  the  high  pressure  results  for 
helium.  Figure  13  shows  the  nitrogen  results  for  a  pressure  of  1  bar  and  frequency  of  8.6  Hz 
which  have  very  clear  similarities  with  the  helium  results  at  5  bar  shown  in  Figure  10. 


8.5  Hz,  1  bar  pressure 


Figure  13,  Fine  Wire  RT  measurements:  Nitrogen,  Ibar,  8,5Hz  for  different  strokes 

(mm) 

With  increasing  pressure  the  nitrogen  results  are  similar  to  the  higher  pressure  helium  results 
with  the  double  minima  feature  being  retained  but  changing  shape  and  showing  a  sharper 
discontinuity.  The  high  frequency  components  and  reduction  in  repeatability  also  set  in,  but 
this  time  increase  to  the  point  where  in  a  few  eases  the  waveforms  becomes  very  unstable  and 
unrepeatable.  For  these  results  the  behaviour  appears  chaotic  and  is  characterized  by  both  the 
presence  of  high  frequency  components  and  a  laek  of  any  repeatable  waveform.  In  the  8.5Hz 
example  shown  in  Figure  14,  there  are  frequeney  eomponents  as  high  as  250Flz.  The  pressure 
variation  for  this  set  of  data  points  is  shown  in  Figure  15.  It  will  be  seen  that  there  is  no  sign 
of  this  instability  in  these  results.  Unfortunately  the  current  limit  imposed  on  the  motor  did 
not  allow  higher  strokes  to  be  investigated  for  this  pressure  and  frequeney. 

It  might  be  expeeted  that  this  unstable  behaviour  might  be  associated  with  some  flow 
transition  and  would  also  be  discernible  in  some  combination  of  stroke  and  pressure  for 
higher  frequeney  data.  This  was  not  found  to  be  the  ease  -  the  only  examples  of  this  extreme 
instability  oecur  in  the  low  frequeney  data  i.e.  at  8.5  or  18  Hz. 

For  the  highest  combinations  of  pressure  frequency  and  stroke  there  is  variability  between 
cyeles,  but  there  are  still  features  of  the  waveform  that  remain  consistent.  Figure  16  shows  the 
waveform  for  nitrogen  for  sueh  eonditions.  It  will  be  seen  that  it  is  remarkably  similar  to  the 
high  pressure,  high  frequency  results  for  helium  shown  in  Figure  12 


8.5  Hz,  18  bar  pressure 


Figure  14,  Fine  Wire  RT  measurements:  Nitrogen,  18bar,  8,5Hz 


8.5  Hz,  18  bar  pressure 
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Figure  15,  Pressure  measurements:  Nitrogen,  18bar,  8,5Hz 

Nitrogen 


Figure  16,  Fine  Wire  RT  measurements:  Nitrogen,  26bar,  43Hz,  stroke  =  14mm 


Summary 

These  results  tend  to  support  the  idea  that  there  are  two  effects  controlling  the  temperature 
behaviour;  Thermal  diffusivity  along  an  axial  temperature  gradient  and  the  transition  from 
laminar  to  turbulent  flow. 

Thermal  Diffusivity 

Firstly  for  both  gases  the  low  stroke  temperatures  appear  to  be  roughly  sinusoidal.  As  the 
stroke  increases  additional  features  begin  to  appear;  initially  in  the  form  of  gradient 
discontinuities,  ultimately  with  the  appearance  in  many  results  of  an  additional  minimum  in 
the  low  temperature  half  of  the  cycle.  The  additional  minimum  always  appears  for  the 
nitrogen  results  but  for  the  helium  results  it  is  absent  in  the  lowest  pressure  measurements. 


A  tentative  conclusion  that  may  be  drawn  is  that  the  thermometer  is  measuring  the 
temperature  of  two  distinct  volumes  of  gas  and  that  the  sudden  changes  in  gradient  that  are 
generally  observable  correspond  to  the  transition  from  one  volume  to  another.  The  interesting 
thing  to  note  is  that  the  transition  becomes  more  abrupt  at  higher  pressures  and  frequencies.  If 
the  mixing  in  the  compression  volume  was  increasing  with  increasing  pressure  and  frequency 
then  it  would  be  expected  that  the  increased  heat  transfer  would  tend  to  smooth  out  the 
transition.  The  observation  that  the  transition  becomes  more  distinct  suggests  that  the  heat 
transport  is  less  effective.  A  possible  explanation  for  this  is  that  the  gas  remains  stratified  and 
that  the  axial  temperature  gradient  is  enhanced  because  the  thermal  diffusion  is  reduced. 
Higher  pressure  reduces  thermal  diffusivity  and  if  the  frequency  is  higher  then  there  is  also 
less  time  for  the  heat  conduction  to  occur.  For  the  low  pressure,  low  frequency,  helium 
results  it  appears  that  the  thermal  diffusion  is  sufficient  to  largely  smooth  out  any  transition. 
This  is  not  seen  in  the  corresponding  nitrogen  results  as  its  thermal  diffusion  is  substantially 
lower  (by  a  factor  of  >8): 

Thermal  Diffusivity  for  Helium  @  1  bar,  300  K  =  1.8  x  10'"^  m^/s 

Thermal  Diffusivity  for  Nitrogen  @  1  bar,  300  K  =  2.1  x  10'^  m^/s 

Transition  from  Laminar  To  Turbulent  Flow 

The  evidence  for  transition  from  a  laminar  to  a  turbulent  regime  is  the  contrast  in  the 
temperature  waveform  stabilities.  In  the  low  pressure  results,  particularly  for  helium,  the 
waveforms  are  very  repeatable  and  show  no  sign  of  randomness  or  instability.  In  the  high 
pressure  results,  particularly  for  nitrogen,  there  is  very  significant  variability  and  in  a  few 
instances  with  high  pressure  nitrogen  we  appear  to  observe  almost  chaotic  behaviour 
displaying  significant  high  frequency  components.  These  results  suggest  that  the  helium 
measurements  are  mainly  in  a  laminar  regime  with  a  tendency  to  turbulence  only  at  the 
highest  pressures.  With  the  nitrogen  results  it  appears  that  they  start  off  in  a  laminar  regime 
but  these  changes  more  quickly  to  being  turbulent. 


7  MEASUREMENTS  MADE  USING  COMPENSATED 
THERMOCOUPLES 

Fast  temperature  measurements  using  a  compensated  thermocouple  technique  have  been 
successfully  implemented  in  IC  engines  (Kar^).  The  technique  has  the  advantage  of  allowing 
substantially  more  robust  sensors  to  be  used  whilst  still  achieving  an  adequate  time  response. 
Although  the  Stirling  cooler  environment  is  less  hostile  than  an  IC  engine,  this  advantage  still 
has  practical  significance  as  the  alternative  of  using  very  fine  wire  sensors  will  always  give 
handling  and  endurance  problems. 

In  the  IC  engine  application  the  temperature  variations  measured  were  generally  quite  large  in 
comparison  with  the  variations  expected  in  the  compression  test  rig  volumes.  As  the  signal  to 
noise  ratio  is  a  significant  factor  in  determining  the  degree  of  compensation  that  can  be 
achieved,  the  challenge  in  this  application  is  to  see  whether  an  adequate  response  is  possible 
with  the  low  signal  levels  expected  in  this  application. 

Description  of  Thermocouples  and  Installation 

In  the  measurement  of  a  fluid  using  a  solid  sensor  there  is  always  a  delay  between  the  change 
in  fluid  temperature  and  its  detection  by  the  sensor.  This  is  due  to  the  time  taken  for  the 
sensor  to  reach  thermal  equilibrium  with  the  gas  and  this  delay  is  determined  by  the  thermal 


mass  and  heat  transfer  rate.  For  fast  response  times  the  sensor  size  is  reduced  to  a  minimum 
so  as  to  improve  the  ratio  of  heat  transfer  to  thermal  mass.  For  sensors  that  are  required  to 
follow  gas  temperatures  with  frequency  components  of  ~  lOOFlz  very  fine  sensors  (such  as  the 
3.8  micron  resistance  thermometer  described  above)  are  required. 

Sensors  that  are  too  large  will  show  an  attenuated  and  delayed  variation  with  higher  frequency 
components  increasingly  filtered  out.  In  principle,  if  the  time  constant  function  is  known  it  is 
possible  to  reconstruct  the  original  gas  temperature  by  compensating  for  the  delay. 

The  compensated  thermocouple  technique  used  here  achieves  this  compensation  by 
positioning  two  thermocouples  with  different  time  constants  such  that  they  see  the  same 
temperature  variation.  The  outputs  from  the  sensors  undergo  signal  processing  to  extract  a 
reconstructed  gas  temperature.  It  was  found  that  the  signal  to  noise  ratio  needed  to  be 
improved  in  order  for  the  technique  to  give  satisfactory  results.  This  was  done  by  sampling 
many  cycles  and  using  ensemble  averaging. 

The  thermocouple  assembly  is  shown  in  Figure  17.  The  two  thermocouples  are  arranged  in 
cross  formation  with  a  separation  of  approximately  0.1  mm.  The  sensor  assembly  is 
positioned  centrally  in  the  top  flange  of  volume  C  such  that  the  ceramic  support  is 
approximately  level  with  the  flange.  The  distance  between  the  thermocouple  junctions  and  the 
top  flange  is  around  3.5  to  4  mm. 

A  comprehensive  report  that  gives  much  more  detail  about  the  probes  and  the  different  signal 
processing  schemes  used  is  attached  as  appendix  2. 


Figure  17,  Pictures  of  thermocouple  probe  (courtesy  of  Prof,  Moriyoshi), 

Left;  3D  view  showing  two  fine  thermocouples  are  in  a  cross  configuration. 

Right:  Side  elevation 


Results  and  Discussion 

Reconstructed  temperatures  were  calculated  using  various  types  of  estimator.  Two  techniques 
using  difference  equations  were  found  to  be  most  satisfactory.  These  were  the  Kalman  filter 
and  the  recursive  least  squares  estimators  and  will  be  denoted  by  ZOH-KF  and  ZOH-LS 
respectively.  Figure  18  shows  the  direct  thermocouple  outputs  with  both  reconstructed 
temperatures  for  one  operating  condition.  The  temperature  variation  for  the  fine  wire 
resistance  thermometer  positioned  at  the  other  end  of  the  regenerator  is  also  shown  for 
comparison. 


Figure  18,  Record  of  the  raw  and  compensated  thermocouple  temperatures  and  cold 
wire  temperatures  (i,e,  at  either  side  of  regenerator).  The  sampling  frequency  was 
decimated  to  1.25  kHz. Helium  gas,  5  bar  fill  pressure,  10  mm  stroke  at  8.6  Hz.  Note  that  the 

two  “R25,  ZOH..”  traces  are  virtually  identical. 


Similar  reconstructed  temperatures  for  a  wide  range  of  operating  conditions  are  given  in 
appendix  2.  It  is  noted  that  for  all  the  results  the  two  reconstructed  temperatures  appear  to 
agree  very  well  with  each  other. 

Although  we  do  not  have  any  means  for  directly  corroborating  the  reconstructed  values  some 
model  tests  were  performed  on  the  reconstructions  and  we  can  also  look  to  see  if  the  general 
features  are  plausible. 

The  model  tests  are  described  in  detail  in  appendix  2.  The  discussion  here  will  be  limited  to 
just  the  physical  implications.  A  few  features  of  the  reconstructed  temperatures  are  worth 
commenting  on; 

•  Firstly  it  is  observed  that  the  phases  of  the  reconstructed  temperatures  are  in  advance 
of  the  raw  thermocouple  data  and  are  close  to  the  phases  of  the  fine  wire  RT 
temperatures. 

•  Secondly  it  is  observed  that  the  peak  to  peak  temperature  variation  is  significantly 
higher  than  for  the  fine  wire  thermometer.  This  suggests  that  the  thermodynamic 
processes  occurring  are  closer  to  adiabatic. 


•  Finally  the  peak  to  peak  temperature  variation  can  be  used  together  with  the  peak  to 
peak  pressure  to  infer  a  mean  polytropic  constant.  These  have  been  calculated  for  a 
number  of  the  helium  results  with  the  finding  that  the  polytropic  constant  is  ~1 .6  for 
low  frequencies  and  pressures  with  the  values  decreasing  as  frequency  and  pressure 
increase.  These  results  suggests  that  there  is  increased  mixing  with  higher  Reynolds 
numbers. 

These  observations  do  not  prove  that  the  technique  has  any  specific  accuracy  but  they  do  at 
least  confirm  that  the  results  are  feasible. 

Conclusions 

The  model  tests  (appendix  2)  and  the  experimental  observations  support  the  general 
conclusion  that  the  temperature  reconstruction  is  useful  provided  that  the  varying  temperature 
is  sufficiently  repeatable  to  allow  ensemble  averaging. 


8  DESCRIPTION  OF  POWER  LOSS  MEASUREMENTS 


The  power  determined  from  the  integral  of  P.dV  in  the  compression  space  is  the  total  power 
input  to  the  gas  and  includes  various  components.  The  purpose  of  this  analysis  is  to  extract 
from  the  power  measurements  the  contribution  that  is  not  accounted  for  in  present  models. 
These  residual  values  can  then  be  investigated  to  see  how  they  vary  with  changing  operating 
parameters. 

The  components  that  make  up  the  total  input  power  to  the  gas  are; 


•  Seal  Loss 

•  Flow  Friction  losses  associated  with  flow  through  the  regenerator 

•  Known  thermal  diffusion  losses  that  correspond  to  boundary  layer  losses  in  the  simple 
gas  spring  configuration.  These  will  be  referred  to  as  Boundary  Gas  Spring  Losses 

•  Residual  values  that  make  up  the  unaccounted  compression  loss.  These  will  be 
referred  to  as  Residual  Compression  Losses 


Seal  Loss  {Psioss) 

Power  is  expended  by  pumping  gas  backwards  and  forwards  through  the  clearance  seal,  ft  can 
be  calculated  from; 


sloss 


dV 

dt 


Where  p  is  the  pressure  difference  and  is  volume  flow  rate  for  the  seal. 

If  the  pressure  variation  is  sinusoidal  with  an  amplitude  pi,  and  laminar  flow  is  assumed  with 
a  constant  radial  clearance  of  T,  then  the  power  expended  is  given  by; 


(Laminar  flow  through  an  annulus  is  assumed) 


P  =  _ E _ 

SIOSS  rs  A  T  -^1 

24//Lp 

where  dp  is  piston  diameter,  Lp  is  piston  length  and  [X  is  gas  viscosity. 

Although  it  is  possible  to  obtain  a  more  accurate  value  by  using  instantaneous  pressure 
values,  this  was  not  considered  necessary  for  these  measurements  and  the  simple  expression 
given  above  was  used  for  all  the  results  given  here. 


Flow  Friction  Losses  {P floss) 

The  power  loss  associated  with  the  flow  through  the  regenerator  is  found  by  evaluating 
f  dV 

P floss  ~  r  P- -  length  of  the  regenerator.  As  the  volume  flow  varies  significantly  from 

^  dt 

one  end  of  the  regenerator  to  the  other,  this  is  done  by  calculating  local  power  loss  values  and 
then  integrating  along  the  length  of  the  regenerator. 

The  power  loss  for  flow  through  an  element  dx  of  the  regenerator  is  given  by 

dP floss  =u{x).pa{x).A^ 

where  pd  is  the  pressure  difference  across  the  element  dx,  u(x)  is  the  flow  velocity  and  is 
the  flow  area. 

For  flow  in  plain  square  mesh  the  friction  factor  generally  has  the  form; 


/  = 


a 


+  b 


where  Nr  is  the  local  Reynolds  number  and  a  and  b  are  constants  denoting  the  laminar  and 
turbulent  components  respectfully 

If  a  linear  velocity  gradient  is  assumed  then  it  is  shown  in  Appendix  1  that  for  sinusoidal 
variations  the  mean  power  loss  is  given  by: 


P  = 

floss 
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Dh-Uj 


d-  !  h  3\  /  A  4\ 

—  .(Mq  +  — .(Mq  ~Ujf) 


where  uo  and  ur  are  the  velocity  amplitudes  at  either  end  of  the  regenerator,  L  is  length  of 
regenerator,  Dh  is  hydraulic  diameter  for  regenerator  mesh  and  0 ,  y/  and  Ud  are  given  by: 


\j/  =  b.p  and  dx  = 


.du(x) 


p  is  gas  viscosity 


The  gas  velocities  at  either  end  of  the  regenerator  can  be  calculated  from  the  compression 
volumes  as  follows: 


For  Volume  C  peak  gas  veloeity  out  of  regenerator  is  given  by  (see  appendix  1) 


1 

{  K  ] 

(dv  ^ 

a 

A 

[K+n+vJ 

y  dt  d 

For  Volume  B  peak  gas  veloeity  into  regenerator  is  given  by 


1 

r  n  ) 

(dV\ 

4' 

\  dt  d 

Va,  Vb,  and  Vc  are  volumes  of  Volumes  A,  B  and  C  respectively.  (Variations  in  pressure  and 
temperature  are  relatively  small.  They  not  affect  volume  flows  very  much  and  were  not  taken 
into  account) 

Finally  the  amplitude  for  the  pressure  difference  across  the  regenerator  can  be  calculated  from 


L 


Dh-Ud 


/  2  2n  ,  /  3  3n 


Values  for  a  and  b  could  have  been  taken  from  the  literature  but  there  is  a  significant  amount 
of  variation  between  different  correlations  and  it  was  decided  to  see  if  values  could  be 
deduced  by  investigating  how  they  affected  the  loss  correlation  obtained. 

It  was  found  that  the  best  correlation  for  nitrogen  gave  a  friction  factor  correlation  of 


/  = 


0.25 


For  helium  the  best  correlation  gave  a  friction  factor  correlation  of 


/  = 


+  b 


where  the  correlation  was  not  sensitive  to  the  value  of  b.  (Note;  The  coefficient  b  relates  to 
the  turbulent  component  of  the  pressure  drop.  For  the  helium  results  the  Reynolds  numbers 
were  low  and  the  pressure  drop  was  dominated  by  the  laminar  component) 

The  correlation  used  to  calculate  the  flow  losses  was 


/  = 


0.25 


This  is  similar  to  expressions  quoted  in  the  literature  for  packed  mesh. 


Boundary  Gas  Spring  Losses  (PBgioss) 

The  losses  assoeiated  with  thermal  diffusion  through  a  boundary  layer  are  the  prineipal  losses 
that  oeeur  in  a  simple  gas  spring  geometry  as  deseribed  in  the  1993  paper  by  Kornhauser  and 
Smith"^  and  investigated  in  Los-FR2.  In  the  more  eomplieated  eonfiguration  investigated  here 
these  losses  will  still  oeeur  wherever  there  is  a  boundary  layer  and  a  pressure  driven  varying 
temperature. 

The  boundary  layer  gas  spring  loss  will  not  be  exaetly  the  same  as  for  the  simple  gas  spring 
eonfiguration  but  it  will  be  assumed  that  the  expression  derived  for  the  gas  spring 
(Kornhauser  and  Smith"^)  will  give  a  reasonable  estimate.  The  equation  used  is: 


'  Bgloss 


(7-1)  1 

^  cosh(y ) .  sinh(y )  -  sin(y  ).  cos(y )  ^ 

Uo9 

7  ' y' 

^  cosh^(y)-sin^(y)  ^ 

Where  po  and  Vo  are  mean  values  of  pressure  and  volume  and  pa  is  the  pressure  amplitude. 
The  parameter  y  is  a  dimensionless  number  that  is  defined  by 


y  = 


^Pe. 


The  Peelet  number  used  here  is  defined  by: 


Pe 


0) 


0^ 

«0 


where  Dh  is  the  hydraulie  diameter:  and  is  the  thermal  diffusivity. 

Using  these  equations  the  Boundary  Gas  Spring  loss  was  ealeulated  for  the  three  separate 
volumes  and  a  total  Boundary  Gas  Spring  loss  was  ealeulated  as  the  sum  of  these 
eomponents. 


PBgioss  (Total)  =  (VolA)  +  (VolB)  +  (VolC) 

It  is  noted  that  the  Peelet  number  for  Volumes  A  and  C  was  typieally  in  the  range  10  to  100 
indicating  that  the  heat  transfer  in  these  volumes  is  tending  toward  the  adiabatic.  In  contrast 
the  Peelet  number  in  Volume  B,  the  regenerator,  is  low  (-0.5)  confirming  that  the  heat 
transfer  regime  is  close  to  isothermal. 


Residual  Power  Loss  (Prioss) 

The  unaccounted  power  loss  is  calculated  from  the  total  power  and  the  components  described 
above: 


TrIoss  Tpn  (Tsloss  T floss  T^Bgloss^ 


Table  3  gives  values  for  a  range  of  operating  parameters  to  give  an  idea  of  the  magnitudes  of 
these  various  components. 


Table  3,  Breakdown  of  Losses 


Operating 

Conditions 

Total  Power 
Input  into 
gas 

Calculated 
Seal  Loss 

Calculated 
Flow  Friction 
Loss 

Calculated 
Boundary 
Gas  Spring 
loss 

Residual 

Loss 

Helium  1  bar 
Stroke  4  mm 
Freq  8.6  Hz 

6.74  mW 

0.075  mW 
1.1% 

2.91  mW 
43.2% 

0.201  mW 
2.98% 

3.56  mW 
52.82  % 

Helium  21  bar 
Stroke  14  mm 
Freq  44  Hz 

3455  mW 

417  mW 
12.07  % 

1088  mW 
31.49% 

13  mW 
0.376% 

1937  mW 
56.06  % 

Nitrogen  1  bar 
Stroke  4  mm 
Freq  8.6  Hz 

6.43  mW 

0.335  mW 
5.21% 

2.74  mW 
42.61% 

0.043  mW 
0.669% 

3.32  mW 
51.63% 

Nitrogen  26  bar 
Stroke  11  mm 
Freq  53  Hz 

4799  mW 

428  mW 
8.93% 

1954  mW 
40.72% 

37  mW 
0.77% 

2379  mW 
49.58% 

Discussion  of  The  Residual  Loss  {P„dLoss) 

The  residual  values  ealeulated  are  aetual  power  losses.  In  order  to  look  at  what  faetors  might 
determine  these  losses  it  is  neeessary  to  transform  them  into  a  non-dimensional  form.  This 
was  done  by  dividing  the  residual  power  loss  by  a  nominal  power  ealeulated  from  the  produet 
of  frequency,  peak  to  peak  pressure  and  swept  volume. 


_  Rloss 


^  ndLoss 


fApAV 


The  values  of  P„dLoss  are  plotted  against  frequency  for  helium  and  nitrogen  at  different 
pressures  in  figures  19  and  20.  These  plots  show  an  obvious  difference  between  the  results 
for  the  two  gases.  For  helium  the  non-dimensional  loss  is  fairly  constant  for  the  1  and  5  bar 
measurements.  For  the  18  bar  measurements  there  is  decrease  in  non-dimensional  loss  with 
frequency.  In  contrast  all  the  nitrogen  measurements  show  an  increase  in  non-dimensional 
loss  with  frequency. 


Helium  Non-Dimensional  Loss  vs  Frequency 
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Figure  19,  Helium:  Variation  of  Non-Dimensional  Loss  with  Frequency  for  Different 

Fill  Pressures 


Ntirogen:  Nondimensional  Loss  vs  Frequency 
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Figure  20,  Nitrogen:  Variation  of  Non-Dimensional  Loss  with  Frequency  for  Different 

Fill  Pressures 


Initial  Investigation  of  Loss  Dependence 

Attempts  to  pursue  a  more  thorough  dimensional  analysis  are  deseribed  below  but  initially  it 
was  found  useful  to  assume  a  dependence  on  Reynolds  number  and  investigate  the  residual 
loss  with  respect  to  the  varied  parameters.  The  Reynolds  number  assumed  is  defined  as; 


where  p  and  p  are  density  and  viscosity  evaluated  for  compression  space,  v  is  piston  velocity 
and  Dh  is  hydraulic  diameter  for  compression  space. 

A  parameter  Z  is  defined  by 


Z  = 


^  ndLoss 


where/  5  and p  are  frequency  stroke  and  pressure  respectively. 

Values  of  Log{z)  were  plotted  against  respective  values  of  Tog(A^).  Indices  kl,  k2  and  k3 
were  varied  to  find  values  that  gave  the  best  correlation  of  the  data  points. 


Helium  Results 

Figure  21  shows  a  graph  of  plotted  against  Log{Nj^)  for  helium  (with  kl,  k2  and 

k3  =  0).  Figure  22  shows  a  graph  of  Log{z)  plotted  against  Log{Nj^)  for  helium.  The  values 
of  kl,  k2  and  k3  that  gave  the  best  fit  are; 

ki=0.3  k2=0.3  k5=0 

I.e.  the  non-dimensional  loss  is  a  function  of  piston  velocity  i.e.  proportional  to  ^ 

It  will  be  seen  that  a  very  respectable  straight  line  fit  is  achieved  from  points  that  originally 
had  a  lot  of  scatter. 


Log(PndLoss)  plotted  against  Log(Nr)  for  Helium 
K1=k2=k3=0 


Figure  21,  Helium:  Log(Z)  plotted  against  Log(Nr)  for  kl=k2=k3=0 


Log(Z)  plotted  against  Log(Nr)  for  Helium 
K1=0.3  K2=0.3  K3=0 


Figure  22  Helium:  Log(Z)  plotted  against  Log(Nr)  for  kl=k2=0,3  and  k3=0 


Nitrogen  Results 

Figure  23  shows  a  graph  of  plotted  against  Log{Nj^)  for  nitrogen  (with  kl,  k2 

and  k3  =  0).  Figure  24  shows  a  graph  of  Log{z)  plotted  against  Log{Nj^)  for  nitrogen  for 
values  of 


kl=  k2  =  0.3  and 
kl=  k2  =  0.5  and 
kl=  k2  =  0.7  and 


k3  =  0 
H  =  0 
H  =  0 


and  the  values  that  gave  the  best  fit: 

kl  =  0.7  k2  =  0.4 


k3  =  0.\ 


The  results  are  plotted  again  for  the  best  fit  in  Figure  25. 

ft  will  be  seen  that  forcing  dependence  on  only  velocity  does  not  give  a  good  correlation.  A 
significantly  better  correlation  is  achieved  with  a  different  dependence  on  the  parameters. 

The  conclusion  to  be  drawn  is  that  the  behaviour  of  nitrogen  and  helium  is  different, 
ft  is  also  noted  that  the  good  correlation  is  not  for  a  single  straight  line,  instead  there  are  two 
regions  which  look  reasonably  linear  but  which  have  different  gradients.  The  Reynolds 
number  at  which  there  is  this  apparent  transition  in  behaviour  is  ~  3000. 

ft  is  noted  that  for  the  nitrogen  results  the  optimum  values  of  kl,  k2  and  k3  were  not  as  clear 
cut  as  for  the  helium  values  but  this  may  well  be  because  we  are  dealing  with  two  regimes 
that  have  a  differing  dependence  on  flow  parameters. 


Log(Z)  plotted  against  Log{Nr)  for  Nitrogen 
K1=0  K2=0  K3=0 


Figure  23,  Nitrogeu:  Log(Z)  plotted  agaiust  Log(Nr)  for  kl=k2=k3=0 


Log(Z) 


Log(PndLoss)  plotted  against  Log(Nr)  for  Nitrogen 
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Figure  24.  Nitrogen:  Log(Z)  plotted  against  Log(Nr)  for  kl=k2=k3=0 
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Figure  25,  Nitrogen:  Log(Z)  plotted  against  Log(Nr)  for  kl=0,7  k2=0,4  and  k3=0,l 


Comments  on  Results 

The  initial  conclusions  to  be  drawn  from  these  results  are; 

•  There  is  a  very  strong  dependence  on  the  Reynolds  number. 

•  Two  regimes  could  be  identified  in  the  nitrogen  results  with  a  fairly  well  defined 
transition  at  a  Reynolds  number  of  3000. 

•  The  helium  results  only  just  extended  up  to  Reynolds  number  values  of  3000  so  it  not 
possible  to  discern  a  transition  at  this  point. 

•  The  dependence  on  frequency  and  stroke  is  quite  marked  for  both  helium  and  nitrogen 
but  the  magnitudes  vary  significantly.  There  appears  to  be  no  dependence  on  pressure 
for  the  helium  results  and  for  nitrogen  the  dependence  was  marginal. 


Nature  of  Residual  Loss 

So  far  we  have  just  referred  to  a  Residual  Loss  without  any  comment  as  to  its  origin.  The 
work  described  does  not  attempt  to  model  a  new  mechanism  but  it  is  nonetheless  desirable  to 
at  least  speculate  as  to  the  nature  of  this  extra  loss. 

In  general  terms  power  can  be  dissipated  in  this  type  of  test  rig  by  three  mechanisms: 

•  Flow  friction 

•  Irreversible  heat  transfer 

•  Work  output 

The  only  significant  flow  friction  mechanisms  are  in  the  seal  and  regenerator  and  these  have 
been  adequately  accounted  for.  The  test  rig  has  only  one  component  that  can  transfer  work 
and  that  is  the  piston  which  is  the  source  of  the  power  input.  It  therefore  looks  likely  that  our 
residual  loss  is  caused  by  an  irreversible  heat  transfer  mechanism  that  is  additional  to  the 
Boundary  Gas  Spring  Loss. 

If  we  look  closely  at  the  tidal  gas  that  spends  some  of  its  time  in  the  compression  space  and 
some  of  the  time  in  the  regenerator  we  see  that  there  is  an  abrupt  change  in  heat  transfer 
conditions  as  gas  moves  from  one  component  to  the  other.  This  discontinuity  will  in  general 
cause  a  significant  irreversible  heat  transfer  and  is  regarded  as  a  strong  candidate  for  our 
missing  loss.  The  mechanism  will  be  briefly  described  here  with  reference  to  the  boundary 
layer  loss  that  is  well  documented. 

Figure  26  shows  a  representation  of  a  simple  gas  spring  at  four  piston  positions.  It  is  assumed 
that  the  boundary  layer  is  relatively  thin  so  that  we  can  regard  the  gas  piston  as  having  a 
lossless  adiabatic  core  with  a  thin  lossy  boundary  layer. 

If  we  start  with  the  piston  midway  in  its  compression  stroke  we  will  assume  that  the  core 
temperature  is  roughly  the  same  as  the  temperature  of  the  walls  and  that  the  heat  transfer  will 
be  close  to  zero  (A).  As  we  proceed  with  the  compression  the  core  gas  temperature  increases 
and  heat  is  transferred  from  the  gas  boundary  layer  into  the  wall.  This  heat  transfer  continues 
during  maximum  compression  Fig  25  (B)  until  we  reach  midway  in  the  expansion  stroke  (C) 


where  the  heat  transfer  is  again  elose  to  zero.  If  we  eontinue  with  the  expansion  then  the  eore 
temperature  beeomes  lower  than  the  wall  temperature  and  heat  transfer  is  reversed  with  heat 
returning  from  the  wall  into  the  boundary  layer.  This  heat  eontinues  until  we  finish  the  cyele 
at  the  midpoint  of  the  eompression  stroke  (A). 

Overall  we  ean  be  sure  of  three  things; 


•  The  mean  gas  temperature  remains  eonstant  henee  there  is  no  net  energy  flow  into  the 
gas. 

•  The  net  heat  transferred  from  the  gas  to  the  wall  is  not  zero.  It  must  be  equal  to  the 
energy  input  into  the  gas  in  the  form  of  work  for  energy  to  be  conserved. 

•  The  power  absorbed  is  a  direct  result  of  the  irreversible  heat  transfer  occurring  in  the 
boundary  layer. 
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Figure  26,  Boundary  Layer  Gas  Spring  Loss  Mechanism, 


If  we  now  look  at  Figure  27  we  see  that  tidal  gas  movement  between  the  compression  space 
(assumed  adiabatic)  and  the  regenerator  is  shown  for  four  points  in  a  cycle.  At  point  A  the 
gas  envelope  is  mainly  contained  within  the  regenerator  and  will  be  at  the  regenerator 
temperature  and  roughly  the  maximum  cycle  pressure.  As  the  pressure  (A  to  C)  decreases  the 
gas  envelope  expands  and  moves  from  the  regenerator  to  the  adiabatic  volume  whilst  also 
tending  to  cool.  While  gas  is  in  the  regenerator  it  will  absorb  heat  and  remain  isothermal;  as 
soon  as  gas  enters  the  adiabatic  volume  it  will  undergo  adiabatic  expansion  and  its 
temperature  will  reduce.  Different  elements  of  gas  will  undergo  differing  amounts  of  cooling 
with  the  result  that  a  temperature  gradient  will  be  set  up  and  irreversible  heat  transfer  will 
occur  with  the  adiabatic  gas  volume.  On  the  compression  stroke  (C  to  A)  the  flow  is  reversed 
and  elements  of  gas  will  be  moving  from  the  adiabatic  volume  into  the  regenerator  whilst 
tending  to  heat  up.  Heat  transfer  within  the  adiabatic  gas  will  continue,  but  in  addition  there 


will  tend  to  be  an  irreversible  transfer  of  heat  between  the  gas  entering  the  regenerator  and 
regenerator  matrix.  Whilst  it  is  possible  that  the  gas  temperature  of  the  adiabatic  gas  is 
synchronised  so  that  it  reaches  the  regenerator  temperature  just  as  it  enters  the  regenerator, 
this  will  not  be  the  general  case. 


Regenerator:  Isothermal  j' 


A. 

B. 

C. 

L' 

/  1 

1 

y~/~/ 

Compression  Space:  Adiabatic 


Figure  27,  Tidal  Gas  Spring  Loss  Mechanism 

We  can  identify  three  distinct  aspects  which  will  tend  to  produce  irreversible  heat  transfer  and 
therefore  compression  losses; 

•  Conduction  of  heat  within  the  adiabatic  volume  driven  by  the  temperature  gradient. 
This  effect  is  analogous  to  the  heat  conduction  in  the  boundary  layer  gas  spring  loss 
and  will  be  governed  by  the  thermal  diffusivity. 

•  Transfer  of  heat  as  above  but  amplified  by  convection.  If  the  gas  remains  stratified  i.e. 
unmixed  then  this  will  not  occur.  As  soon  as  there  is  significant  mixing  then  this  effect 
is  likely  to  become  dominant  in  the  way  that  turbulent  heat  transfer  is  dominant  in 
steady  flow  heat  transfer. 

•  Transfer  of  heat  from  gas  to  regenerator  matrix  as  gas  enters  the  regenerator. 

From  the  above  it  is  seen  that  there  is  an  obvious  mechanism  for  irreversible  heat  transfer  and 
with  two  additional  components,  it  is  substantially  more  complicated  than  the  mechanism 
occurring  in  a  simple  Boundary  Layer  gas  spring.  An  appropriate  term  for  this  effect  is  a 
'^Tidal  Gas  Spring  Loss'" 


Dimensional  Analysis 


The  correlations  obtained  above  are  surprisingly  good  if  they  really  represent  the  loss  outlined 
above.  However,  to  be  of  any  real  use  it  is  desirable  to  formulate  these  correlations  in  a  non- 
dimensional  form.  As  a  starting  point  an  analysis  is  pursued  with  the  more  obvious 
parameters  listed  below  in  table  4 

Table  4,  Parameters  used  in  Dimensional  Analysis 


Variable 

Symbol 

Units 

Frequency 

f 

s-' 

Viscosity 

fi 

kg.m''.  s’' 

Density 

P 

kg.m'^ 

Hydraulic  Diameter 

Dh 

m 

Velocity 

V 

m.  s'* 

Thermal  Diffusivity 

a 

m**.  s'* 

Gamma  (ratio  of  specific  heats) 

7 

none 

The  ratio  of  specific  heats  is  dimensionless  and  is  assumed  to  form  its  own  non-dimensional 
group.  There  are  three  fundamental  dimensions;  mass  (kg),  length  (m)  and  time  (s).  There  are 
assumed  to  be  six  variables,  hence  there  will  be  three  non-dimensional  groups. 

Repeating  variables  are  chosen  to  be  density  p,  velocity  v  and  hydraulic  diameter  Dh.  The 
non-dimensional  groups  P{1),  P{2)  and  P{3)  will  then  have  the  forms; 


p(l)  = 


P{2)  =  p  .V  .D,  .a 


P{3)  =  p‘‘\v'’\Df 


The  indices  are  found  by  equating  dimensions,  the  three  dimensionless  groups  are; 


P{\)  = 


P-V-Dh 

P 


P{2)  = 


P{3)  = 


f-Du 


Hence  the  non-dimensional  loss  factor  is  expected  to  have  a  form; 


Pn,Loss  =  K.f{r)\ 


P-V-Dh 

V  P  ) 


\  a  ) 


V  V  ; 


+  other  terms 


These  dimensionless  groups  can  be  identified  as; 


p(l)^AV.A 

Reynolds  Number 

P 

P(2)  = 

Peclet  Number 

A7, 

a 

Strouhal  Number 

V 


The  formulation  of  Peclet  number  ean  be  improved  to  avoid  disguising  the  density 
dependenee  in  the  thermal  diffusivity.  The  thermal  diffusivity  is  defined  by: 

.  .  ,  K 

a  IS  given  by  a  = - 

P-Cp 

K  is  thermal  eonduetivity  and  Cp  is  speeifie  heat  eapaeity 
If  we  substitute  this  in  the  expression  above  we  obtain 

^(2)  =  ^-^.p.Cp  =  p.v.D^.^ 

Comparing  this  with  the  Reynolds  number  we  see  that  both  dimensionless  groups  have  the 
variables  p,  v  and  Dh  to  the  same  power  and  that  the  other  parameters  are  substantially 
eonstant  for  the  gases  eoneerned.  The  result  is  that  for  individual  gases  the  Peclet  number  and 
Reynolds  number  will  exhibit  the  same  parameter  dependence  -  in  effect  we  only  have  one 
dimensionless  group.  In  fact  the  ratio  of  Reynolds  number  to  Peclet  number  defines  another 
dimensionless  group  that  is  frequently  used  in  heat  transfer  problems,  the  Prandtl  number  Np/. 

N  =N  *N 

The  Prandtl  is  typically  equal  to  0.7  for  most  gases.  The  difference  between  the  values  for 
nitrogen  and  helium  is  small  and  would  not  account  for  significant  differences  in  their 
behaviour. 

If  we  now  attempt  to  use  these  dimensionless  groups  to  correlate  the  helium  data  we  hit  an 
immediate  problem.  In  the  initial  analysis  given  above  we  found  that  the  non-dimensional  loss 
was  proportional  to  the  Reynolds  number  and  The  Reynolds  number  is  unchanged  and 
we  need  to  find  a  combination  of  the  dimensionless  groups  that  will  give  the  required  velocity 
dependence.  With  the  only  the  three  groups  given  here  this  does  not  seem  possible  and  the 
conclusion  to  be  drawn  is  that  we  are  missing  at  least  one  important  variable. 

If  the  Mach  number  is  invoked  this  will  allow  the  velocity  dependence  to  be  expressed  in 
terms  of  the  velocity  of  sound.  However  for  the  velocities  encountered  in  this  investigation 
compressibility  effects  appear  fairly  unlikely. 

In  the  discussion  above  concerning  possible  loss  mechanisms  in  tidal  gas  flow  we  did  not 
consider  how  the  phase  between  pressure  variation  and  mass  flow  might  influence  the  loss. 
The  likelihood  is  that  it  is  a  significant  parameter  and  it  may  be  that  the  apparent  dependence 
on  velocity  is  an  expression  of  this.  Taking  the  modelling  or  dimensional  analysis  further  to 
investigate  this  is  beyond  the  scope  of  the  work  reported  here  and  we  will  have  to  conclude 
that  whilst  there  does  seem  to  be  effects  that  can  be  well  correlated  the  determining 
parameters  have  yet  to  be  fully  identified. 


9  OVERALL  SUMMARY  AND  CONCLUSIONS 


The  most  interesting  feature  of  this  work  is  the  general  agreement  between  the  resistanee  wire 
thermometry  results  and  the  power  loss  measurements.  Both  sets  of  results  exhibit  ehanges  in 
flow  behaviour  eonsistent  with  the  idea  of  a  laminar  to  turbulent  transition.  Furthermore  a 
region  of  transition  defined  by  a  Reynolds  number  suggested  by  the  power  measurements  is 
reasonably  supported  by  the  thermometer  results.  In  addition  the  resistanee  wire  thermometry 
results  ehanges  in  temperature  profile  are  eonsistent  with  thermal  diffusivity  effeets.  It  is  also 
elear  that  the  gas  is  quite  stratified  in  the  low  pressure  results  and  that  this  stratifieation  is  still 
present  to  a  lesser  extent  in  the  high  pressure  results  that  are  thought  to  turbulent. 

A  simple  dimensional  analysis  eonfirms  that  the  thermal  diffusion  effeets  should  be 
determined  by  a  Peelet  number  whilst  the  transition  from  laminar  to  turbulent  should  be 
determined  by  a  Reynolds  number.  For  the  Peelet  and  Reynolds  numbers  defined  it  is  found 
that  they  have  the  same  dependenee  on  the  parameters  varied  i.e.  pressure,  frequeney  and 
stroke.  The  use  of  gases  with  different  physieal  properties  should  allow  differentiation  of 
these  groups  but  this  could  not  be  done  here  with  the  limited  results  available 

The  more  detailed  conclusions  that  can  be  drawn  are; 

•  Cyclic  temperature  measurements  were  made  for  the  gas  in  the  fixed  cylindrical  volume 
(Vol  C)  attached  to  the  end  of  the  regenerator.  The  thermometer  used  was  of  a  novel  type 
that  used  two  relatively  slow  thermocouples  and  a  compensation  technique  to  reconstruct 
the  actual  gas  temperature.  The  measurements  made  with  the  compensated  thermocouples 
did  not  give  any  information  about  the  compression  losses  occurring  in  this  gas  spring 
configuration.  However  they  did  show  that  this  technique  could  be  considered  in  this  type 
of  application,  provided  that  the  noise  levels  were  low  enough.  For  the  measurements 
described  here  the  noise  could  only  be  sufficiently  reduced  by  ensemble  averaging.  This 
would  be  a  limitation  in  circumstances  where  the  waveform  is  not  very  repeatable. 

•  Cyclic  temperature  measurements  were  made  for  the  gas  in  the  top  of  the  compression 
volume  (Vol  A)  adjacent  to  the  end  of  the  regenerator.  The  thermometer  used  was  a 
resistance  thermometer  that  used  very  fine  wire  to  give  a  rapid  response. 

•  The  temperature  profiles  measured  showed  a  range  of  shape  and  repeatability  depending 
on  the  gas  used  and  also  on  the  stroke,  frequency  and  pressure. 

•  A  significant  part  of  the  stroke  dependence  appears  to  be  due  to  the  portion  of  time  in 
which  the  thermometer  was  sensing  the  tidal  gas  flow  and  the  proportion  attributable  to 
the  gas  which  always  remains  in  the  compression  space. 

•  This  sensing  of  two  distinct  gas  volumes  is  supported  by  the  apparent  gradient 
discontinuities  that  can  be  seen  in  the  profiles.  There  are  no  corresponding  discontinuities 
in  the  pressure  waveforms. 

•  The  variation  in  the  gas  temperature  waveforms  appear  to  be  due  to  two  effects: 

o  Changes  in  the  axial  temperature  gradient  as  the  thermal  diffusion  varies, 
o  A  tendency  with  increased  pressure  towards  waveform  variability,  instability  and 
the  presence  of  high  frequency  components,  this  suggests  a  transition  to  a  more 
disordered  regime:  perhaps  Laminar  to  Turbulent. 


•  The  helium  results  tended  to  mainly  oeeupy  the  ordered  laminar  regime.  The  nitrogen 
results  appeared  to  straddle  both. 

•  A  few  high  pressure,  low  frequeney  nitrogen  results  showed  extreme  instability  that  eould 
be  deseribed  as  ehaotie  -  no  repeatable  waveform  was  diseernible.  It  is  possible  that  these 
results  just  happened  to  straddle  a  laminar/turbulent  transition  but  with  the  limited  data 
available  this  eonelusion  must  be  regarded  as  tentative. 

•  The  power  loss  measurements  showed  that  a  loss  existed  that  eould  be  eorrelated  by  a 
eombination  of  Reynolds  number  and  the  experimental  parameters  pressure,  frequeney 
and  stroke.  Attempts  to  formulate  this  in  a  non-dimensional  form  were  not  sueeessful 
suggesting  that  key  parameters  may  not  have  been  invoked.  One  effeet  in  partieular  that 
needs  eloser  attention  is  the  effeet  of  phase  angle  between  mass  flow  and  pressure 
variation. 


10  FUTURE  WORK 

For  this  investigation  into  eompression  losses  to  be  eontinued  two  areas  of  work  stand  out: 

•  Further  experimental  measurements  repeating  resistanee  wire  thermometry  and  power  loss 
measurements  for  different  gases  and  operating  eonditions.  These  would  eoneentrate  on 
the  transitional  regimes  already  identified  to  provide  a  better  experimental  deseription. 

•  A  goal  to  establish  a  model  for  these  proposed  Tidal  gas  Spring  Losses  appears  to  be  very 
important  aspiration.  It  will  remain  diffieult  to  make  mueh  sense  of  the  results  without  a 
model  to  aet  as  a  framework. 
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APPENDIX  1:  PRESSURE  DROP  IN  REGENERATOR 


BACKGROUND  EQUATIONS 
Kays  and  London 

Shear  Stress  given  by: 


r  = 


l.p 


/is  friction  factor  G  is  mass  flow  rate  per  area  (i.e  mass  flux)  =  u.p 
p  is  fluid  density 
G  =  u.p  where  u  is  fluid  velocity 

f.p.u^ 


dP.A^  =  t.A^, 


Ax  is  flow  area 


r  = 


A.  = 


2 

Vol 

L 


A...  = 


4.V0I 


A.  D, 


A„  is  wetted  area 


Vol  is  wetted  volume 


Dh  is  hydraulic  diameter  L  is  regenerator  length 

f.p.u^  4.L 


dP  =  T.^  =  T.— 
A.  D, 


dP  = 


dP  = 


{f.p.u^).2.L 


For  flow  through  mesh  friction  factor  has  form 


/  = 


+  b 


N,= 


P-U-Dh 

P 


APPLICATION  TO  COMPRESSION  LOSS  RIG 
Power  Loss 

Flow  velocity  and  hence  pressure  drop  varies  along  length  of  regenerator.  It  is  therefore 
necessary  to  evaluate  pressure  drop  power  loss  for  small  length  dx  and  then  integrate  over 
length  of  regenerator.  It  will  be  assumed  that  flow  velocity  varies  linear  along  length  of 
regenerator 

X 

u(x)  =  Uq - .Uj  for  x=L  u(x)=ul 

L 

uq  is  velocity  at  entrance  to  regenerator,  Ud  is  velocity  difference  across  regenerator 

For  length  of  regenerator  dx  with  fluid  velocity  u(x)  and  pressure  drop  Pd(x),  power  dissipated 

is  given  by: 


dW  =  u(x).P^(x).A^ 


Substituting  expressions  from  above 

Pdi^)  =  i 


p.u{x).D^ 


J 


p.u{x)  .l.dx  _  2 


a.p.u{x) 

~~D,, 


+  b.p.u{xf 


.dx 


dW  = 


2.A^ 


a.ju.u{xy 

V  D, 


+  b.p.u{xy 


.dx 


Let  6  =  ,  y/  =  b.p  and  dx  =  — —.du{x) 


Then  Power  integrated  over  length  of  regenerator  is  given  by: 


I  [d.M(xy  +  y/.M(xy  )du{x)  ■■ 


2.A^  L 


^d 


6.u{xy  ^  y/.u{xy 


evaluated  over  veloeity  limits  uo  and  ul 


W  = 

peak 


Dj^Md  L 


y-(Wo  -ul)  +  ^.{ul  -ul) 


Where  uo  and  ul  are  the  peak  veloeity  amplitudes  at  either  end  of  the  regenerator 
The  above  doesn’t  take  aeeount  of  the  time  variation.  In  prineiple  this  expression  ean  be 
integrated  for  sinusoidal  veloeities.  However  as  the  pressure  drop  inereases  the  veloeity  will 
reduee  henee  it  is  eonsidered  that  the  power  variation  will  be  rounded  off  and  quadratie 
dependenee  on  power  is  a  good  approximation.  The  mean  power  is  then  related  to  the 
maximum  instantaneous  power  by: 

Mean  power  =  0.5  *(maximum  instantaneous  power) 

Henee 


=  ■ 


A-L 

Dh 


/  3  3n  ,  /  4  4n 

—  ■{Uq -Ul)  +  —-{Uq  -Ul) 


Estimation  of  Peak  Gas  Velocities 

Total  of  gas  is  assumed  to  be  eonstant  over  eyele  i.e.  leakage  past  seal  is  negleeted 

M.P.V 

m^ot  =  m  = -  for  eaeh  gas  volume 

R.T 

M  is  mass  of  1  mole  P  is  pressure  V  is  volume 
R  is  universal  gas  eonstant  T  is  temperature 

We  need  to  evaluate  the  peak  values  of  max  volume  flux  V  at  the  boundaries  between  the 
volumes.  The  peak  veloeities  ean  then  be  ealeulated  from 


The  volume  fluxes  are  ealeulated  from  the  mass  variations -  in  eaeh  volume. 

dt 


pr  _  i  dm 
p  dt 

For  the  volumes  that  are  eonstant 


dt  R.T  dt 


V-—  ^ 

p  R.T  dt  P  dt 


Effeet  of  temperature  variation  is  assumed  to  be  small  and  mean  temperature  for  all  volumes 
at  maximum  gas  veloeity  is  assumed  to  be  equal. 

Pressure  variation  is  driven  by  volume  variation  eaused  by  piston  ean  be  ealeulated  from 


\K1 


RT)\V^+V,+V^ 


Differentiating  the  only  variable  is  eompression  volume  Va 


^  =-p. 


K+K+K 


For  Volume  e  peak  gas  veloeity  out  of  regenerator  is  given  by 


A  [K+K+KA  dt 


For  Volume  b  peak  gas  veloeity  into  regenerator  is  given  by 


A  dt 


Pressure  Drop  across  Regenerator 

Total  Pressure  drop  is  given  by  integral  of  pressure  drop  along  length  of  regenerator 
Rd (^)  “  {6.u{x)  +  \f/.u{xY^dx 

D  h 


AP  =  — 


Dh  .u^ 


•j  {o.uix)  +  if/.u{xY  )du(x) 


2.L  6.u{xY  y/.u{xy 

Dh-Ud  2  3 


evaluated  over  veloeity  limits  m  and  ul 


-(^0  ^l)'^  ~  -(^0  ^l) 

Dh-Ud  L2  3 


R  \d-  .  ,2  , ,2^  ,  ¥ 


APPENDIX  2:  REPORT  ON  COMPENSATED  THERMOCOUPLE 
RESULTS 


1  Thermocouple  compensation 

1,1  Introduction 

Two  fine  wire  thermocouples  of  different  diameters  were  installed  after  the  regenerator  exit. 
These  are  R  type  thermocouples  (Platinum/Pt-13%  Rhodium)  with  nominal  diameters  of  25 
and  50  microns.  The  junctions  were  specially  made  by  percussion  welding.  Traditionally,  the 
junctions  were  created  by  spot  welding,  which  has  a  diameter  2-2.5  times  bigger  than  the 
diameter.  The  bigger  the  diameter  of  a  junction,  the  bigger  the  thermal  inertia,  and  so  reduces 
the  time  response.  Percussion  welding  can  form  a  junction  no  bigger  than  the  original 
diameter,  realizing  the  full  potential  of  the  fine  wires.  Each  welded  thermocouple  was 
inspected  under  the  microscope.  The  wire  diameter  has  to  be  uniform  within  10%  for  it  to  be 
accepted. 

The  fine  wires  were  supported  by  four  prongs  made  of  the  same  thermocouple  materials  as 
shown  in  Figure  1.  The  wires  are  arranged  in  a  X  configuration  with  the  junctions  situated  at 
the  centre  of  the  probe.  The  finer  wire  is  raised  slight  above  the  thicker  wire,  so  that  the 
junctions  are  separated  by  0.1  mm. 


Figure  1  Pictures  of  thermocouple  probe  (courtesy  of  Prof.  Moriyoshi).  Left:  3D  view  showiug  two  flue 
thermocouples  are  iu  a  cross  coufiguratiou.  Right:  Side  elevatiou 

Even  though  very  fine  thermocouples  were  used  in  the  experiment,  the  time  response  was  still 
insufficient  to  measure  the  gas  temperature  accurately.  Figure  2  illustrates  this  point.  The 
temperature  variations  in  a  cycle  are  plotted.  TC25  and  TC50  represented  the  25  and  50 
micron  thermocouple  temperatures  respectively.  The  cold  wire  temperature  is  also  plotted  for 
comparison.  The  cold  wire  was  installed  in  the  compression  space,  4.5  mm  upstream  of  the 
regenerator  entry.  Since  the  response  of  the  cold  wire  is  at  least  100  times  faster  than  the 
thermocouple,  the  measured  temperature  may  be  regarded  as  instantaneous. 

The  beginning  of  the  cycle  is  defined  as  zero  piston  displacement  when  the  piston  is  moving 
up.  The  cold  wire  temperature  peaks  at  about  1/4  of  the  cycle,  when  the  piston  has  the 
maximum  displacement.  This  means  the  temperature  rise  is  mainly  due  to  compression. 


However,  the  peak  temperatures  of  TC25  and  TC50  occur  later  and  at  different  parts  of  a 
cycle  with  different  magnitude.  TC50  has  smaller  peak  and  lags  behind  TC25.  The  results 
strongly  indicate  that  the  TC25  and  TC50  are  not  fast  enough  to  response  to  the  changes  in 
the  gas  temperature.  This  section  will  explain  the  theory  of  thermocouple  compensation 
techniques  and  how  they  are  being  applied  to  the  measurements.  The  quality  of  compensation 
will  be  evaluated  by  model  validation. 


No.  of  cycles 

Figure  2  Typical  results.  Helium  gas,  5  bar  fill  pressure,  10  mm  stroke  at  36  Hz. 


1,2  Theory 

A  well  designed  thermocouple  (convective  heat  transfer  is  dominant  compared  to  conductive 
and  radiative  heat  transfer)  can  be  modelled  as  a  first-order  lag  system,  in  which  the  time 
constant  (r)  defines  how  fast  a  thermocouple  will  response  to  a  change  in  temperature.  The 
response  may  be  described  mathematically  as _ 


^  ^  dr 

7;=r+r— 

^  dt 

1) 

where; 

Tg  is  the  temperature  of  the  environment 
T  is  the  thermocouple  measured  temperature 
t  is  time 

Hoshino  [1]  investigated  the  effect  of  conduction  on  the  thermocouple  temperature  under 
Stirling  engine  operating  conditions.  The  simulation  results  have  shown  that  if  the  length-to- 
diameter  ratio  {Lid)  of  a  thermocouple  is  bigger  than  40-80,  the  conduction  effect  is 
negligible.  The  exact  value  depends  on  the  convective  heat  transfer  rate.  The  lowest  Lid  of  the 
thermocouple  probe  is  65,  so  the  conduction  effect  is  expected  to  be  small.  Given  the 
maximum  thermocouple  temperature  swing  observed  is  about  10  K,  and  the  regenerator 


housing  and  other  connected  components  were  more  or  less  at  the  average  thermocouple 
temperature.  There  was  only  ~5  K  temperature  difference  at  about  300  K  to  drive  the  radiative 
heat  transfer  process.  Therefore,  the  radiation  heat  transfer  is  also  expected  to  be  small. 

When  the  gas  temperature  changes  with  time,  Eq.  1  states  that  the  measured  and  the  gas 
temperatures  will  be  different.  This  difference  is  called  the  dynamic  error  due  to  finite  heat 
transfer  between  the  gas  and  the  thermocouple.  The  dynamic  error  can  be  minimised  by 
having  smaller  thermocouples,  and  materials  with  smaller  density  and  specific  heat,  but  it 
cannot  be  totally  eliminated.  In  the  process  of  thermocouple  compensation,  the  time  constant 
is  not  minimised,  instead  it  is  being  estimated  by  some  means.  Then  the  gas  temperature  can 
be  back  calculated  from  the  measured  temperature  and  its  time  derivatives.  Unfortunately,  the 
time  constant  is  not  a  property  of  the  thermocouple,  but  depends  on  parameters  that  can  affect 
the  heat  transfer,  such  as  gas  density,  Reynolds  number  and  thermal  conductivity,  etc. 

The  time  constant  may  be  estimated  analytically  using  established  Nusselt  correlations  for  the 
appropriate  geometry,  e.g.  Collis-William’s  equation  [2].  This  method  has  limited  use 
because: 

•  Nusselt  correlations  are  derived  in  steady  flow,  so  they  may  not  work  in  an  oscillatory 
flow  as  found  in  Stirling  engines. 

•  Flow  velocity  is  required,  which  is  often  not  available. 

•  Some  fluid  and  material  properties  are  functions  of  temperature. 

The  alternative  is  to  measure  the  time  constant  by  creating  a  step  change  in  thermocouple 
temperature  in  the  test  flow.  By  observing  the  temperature  decay,  the  time  constant  is  the  time 
required  for  the  junction  temperature  to  decease  to  63.2%  of  the  difference  between  the  initial 
and  final  temperatures.  The  step  change  is  usually  achieved  by  passing  an  electric  current 
through  the  thermocouple  wire.  Some  researchers  [3]  found  that  the  heat  transfer 
characteristic  is  different  when  the  thermocouple  is  subject  to  internal  heating  rather  than 
external  heating  (when  it  is  working  as  a  thermometer).  Since  the  time  constant  varies  with 
flow  conditions,  it  needs  to  be  measured  for  every  test  flow  condition  and  this  can  be  rather 
time  consuming.  When  the  flow  is  constantly  changing  with  time,  only  an  average  time 
constant  can  be  measured.  The  accuracy  of  the  compensation  will  be  determined  by  the  range 
of  the  time  constant.  It  may  be  acceptable  for  some  applications,  but  may  be  not  for  others. 

The  latest  development  has  been  in  two-thermocouple  techniques.  These  techniques  allow  the 
time  constant  to  be  estimated  in  situ.  No  heating  or  cooling  is  involved.  In  some  techniques 
the  instantaneous  time  constants  can  be  estimated,  as  opposed  to  an  average  value  from  the 
empirical  method.  The  key  assumption  of  these  techniques  is  that  the  two  thermocouples  are 
subjected  to  the  same  temperature.  This  may  be  violated  in  the  case  of  a  steep  temperature 
gradient.  The  model  validation  section  later  on  will  verify  this  assumption.  Many  two- 
thermocouple  techniques  have  been  proposed:  Cambray  [4],  Tagawa  et  al.  [5,  6,  7],  Forney 
and  Fralick  [8],  Hung  et  al.  [9,  10]  and  Kar  et  al.  [11].  This  report  will  only  give  a  brief 
overview  of  these  techniques.  For  more  details,  please  refer  to  [1 1,  12  and  13].  The  most 
important  distinction  between  these  techniques  is  whether  the  time  constant  is  assumed  to  be 
time  invariant.  A  time  invariant  assumption  requires  the  environment  conditions  not  to  change 
with  time. 

Fet  subscript  1  and  2  denote  the  first  and  second  thermocouples  respectively.  The 
thermocouple  equations  are: 


There  is  a  singularity  in  this  solution  when  the  denominator  approaehes  zero.  Moreover,  the 
evaluation  of  the  time  derivatives  will  invariably  amplify  any  measurement  noise  in  the 
signals,  this  will  tend  to  make  the  equation  ill-defined  over  a  region.  The  result  is  unphysieal 
and  erratie  time  eonstant  estimates.  The  measurement  noise  ean  be  handled  by  least  square 
methods  [5,  6,  9,  13].  Instead  of  equating  Eq.  2  and  3,  these  methods  seek  to  minimise  the 


differenee  between  and  T„2  in  a  least  square  sense  in  N  number  of  samples 


This  oeeurs  when 


For  whieh  the  solution  for  the  time  eonstants  is 


In  Eq.  8,  the  following  abbreviation  has  been  used  to  make  the  equations  easier  to  interpret: 


1 

N 


;  Ar  =  r,-7;;  G, 


/=! 


— -  and  G, 
dt  " 


dt 


To  traek  the  ehanges  in  the  time  eonstants,  N  is  set  to  be  smaller  than  the  total  number  of 
samples.  It  is  apparent  that  the  ehoiee  of  is  a  eompromise  between  fast  tracking  and  noise 
rejection.  In  the  limits  when  A^=  1,  the  estimation  become  a  recursive  algorithm.  O’Reilly  et 
al.  [14],  Kar  et  al.  [11],  and  Grys  and  Minkina  [15]  applied  various  recursive  algorithms  for 
time  constant  estimation.  The  major  difference  in  the  algorithms  lies  in  which  estimators  they 
used.  For  example,  Kar  et  al.  used  a  Kalman  filter,  while  Grys  and  Minkina  employed  the 
popular  least  mean  squares  and  recursive  least  square  estimators.  Since  recursive  algorithm 
will  be  applied  in  this  work,  their  theories  will  be  covered  fully  in  Analysis  method  section. 


1,3  Scope 

This  work  focuses  on  both  nitrogen  and  helium  data  at  a  fixed  stroke  of  10  mm.  Three  filling 
pressures  (1,  5  and  18  bar)  and  three  frequencies  (8.6,  18  and  36  Hz)  have  been  investigated, 
that  is  a  total  of  9  operating  points  for  eaeh  gas.  A  longer  stroke  is  not  eonsidered,  as  there  is 
no  experimental  data  for  some  pressures  and/or  frequeneies  due  to  limits  on  the  maximum 
input  voltage  or  current  to  the  compressor. 


1,4  Analysis  Methods 

The  general  flow  of  the  analysis  method  is  illustrated  in  Figure  4.  The  amplified  thermocouple 
voltages  were  digitised  and  reeorded  using  the  Flyde  data  aequisition  system.  As  deseribed  in 
the  data  proeessing  seetion,  the  voltages  are  eorreeted  for  amplifier’s  zero  voltages  and 
sampling  delay,  and  finally  ealibrated  to  temperatures.  The  ealibrated  temperature  has  an 
uncertainty  of  0.2  K. 


1,4,1  Ensemble  averaging 

The  thermoeouple  has  a  low  sensitivity  (about  5  pV/K).  To  improve  the  signal-to-noise  ratio, 
ensemble  averaging  has  been  employed.  It  has  been  found  that  interference  espeeially  from 
the  mains  is  effectively  eliminated.  However,  the  results  will  be  limited  to  “average” 
eonditions;  cyelie  variability  cannot  be  investigated.  This  would  not  be  a  serious  limitation  as 
the  measured  temperatures  are  highly  repeatable,  e.g.  see  Figure  3. 


Figure  3  Thermocouple  temperatures  (TC25)  for  360  cycles,  5  bar  fill  pressure,  10  mm  stroke  at  36  Hz. 

The  ensemble  averaged  temperatures  are  extended  by  joining  the  end  of  cycle  to  the 
beginning  three  times,  so  that  there  are  3  cyeles  of  data.  There  are  several  reasons  for  doing 
so; 


•  Some  filters  have  lead-in  and  lead-out  periods,  in  whieh  the  filtered  signal  is 
meaningless.  Extending  the  data  will  give  all  the  lead-in  in  the  first  eyele  and  the  lead 
out  in  the  third  eyeles,  leaving  the  seeond  eyele  free  of  filtering  errors. 

•  At  high  frequeney  oyelie  frequeney  eonditions,  there  are  not  enough  data  points  in  a 
eyele  for  the  Gaussian  smoothing  filter  to  work. 

•  The  time  derivatives  of  temperatures  are  usually  ealeulated  by  eentral  differenee  (3- 
point  seheme).  This  is  not  possible  for  the  first  and  the  last  data  points.  Having  two 
redundant  eyeles  will  overeome  this  problem. 

A  fixed  number  of  samples  were  eolleeted  for  eaeh  experiment,  irrespeetive  of  the  eyelie 
frequeney.  Therefore,  the  higher  the  eyelie  frequeney,  the  fewer  data  points  eaeh  eyele  has, 
and  the  greater  the  number  of  eyeles.  So  at  the  lowest  frequeney  (8.6  Hz),  there  are  only  85 
eyeles.  At  36  Hz,  the  number  of  eyeles  inereases  to  360.  This  effeet  will  be  diseussed  in  detail 
later. 


1,4,2  Time  derivatives  and  low  pass  filtering 


The  time  derivatives  of  temperatures  are  needed  for  some  estimators  and  for  temperature 
reeonstruetion.  If  low  pass  filtering  is  not  required,  they  are  evaluated  by  eentral  differenee. 
When  low  pass  filtering  is  desired,  the  temperature  signals  are  filtered  by  a  Gaussian 
smoothing  filter.  The  filter  suppresses  the  high  frequeney  noise  by  eonvolving  the  signal  with 
a  Gaussian  kernel.  The  Gaussian  kernel  is  defined  by; 


g{t)  =  exp[- 


9) 


where  /I  is  a  smoothing  parameter  that  sets  the  eut-off  frequeney.  Let  the  thermoeouple  signal 


be  x(t) 

;  the  fdtered  signal  y(t)  is  then: 

yit)  =  x(t)  0  git) 

10) 

The  time  derivatives  ean  be  found  easily  by  using  the  eommutative  law  of  eonvolution  as 
follows: 


dt  dt 

11) 

dt 

where  dg(t)/dt  =  -2/ltexp(-/lt^). 


Compensated 

temperatures 


Figure  4  Analysis  method  flow  chart 


1,4,3  Estimators  and  model  structure 

It  is  unclear  if  the  ehange  in  pressure  and  veloeity  in  the  regenerator  exit  ean  eause  strong 
variations  in  the  time  eonstants.  If  so,  the  two  thermoeouple  teehniques  must  be  able  to  traek 
the  changes  by  means  of  reeursive  algorithms.  However,  sinee  these  teehniques  have  to  plaee 
more  emphasis  on  more  reeent  samples  than  past  samples,  they  are  inherently  more  prone  to 
noise  and  have  bigger  error  (bias)  in  the  estimation.  Thus,  in  some  ease,  an  average  time 
eonstant  with  small  bias  ean  be  more  aeeurate  than  a  varying  time  eonstant  with  a  big  bias.  It 
is  therefore  proposed  to  estimate  an  average  time  eonstant  first.  This  will  establish  the 
benehmark  with  whieh  the  time  varying  results  ean  be  eompared.  In  this  way,  any 
improvement/deterioration  in  aecuraey  ean  be  assessed.  Only  if  the  time  varying  results  show 
improvement  in  the  model  validation  will  they  be  aeeepted. 


The  response  equations  (Eqs.  2  and  3)  ean  be  expressed  in  a  matrix-veetor  form: 


m)-T,{k)= 

1 - 1 

1 

1 _ _ _ 1 

h\ 

A2J 

12) 

y  = 

x^e 

where  k  is  the  k**'  sample  in  a  diseretised  signal.  The  left  hand  side  of  the  equation  is  ealled 
system  output  (y).  The  right  hand  side  eonsists  of  a  regressor  (x)  and  parameter  veetor  (0). 
This  general  form  is  well  known  in  the  system  identifieation  literature  [16,  17].  Onee  the 
thermoeouple  equations  are  in  this  form,  the  time  eonstants  ean  be  estimated  by  different 
estimators,  sueh  as  a  Kalman  filter.  Eq.  12  demonstrates  one  model  strueture  in  whieh  the 
differential  equations  are  used  direetly.  Alternatively,  the  thermoeouple  equations  ean  be 


model 

led  by  differenee  equations  [9],  whieh  yield  another  model  strueture: 

T,{k)-T,{k-\)  =  ^ 

-  T,{k)-T,{k-\)  ~ 
j,{k-\)-m-\)_ 

T 

l-exp(-r7ri) 

13) 

where  4  is  the  sampling  time  whieh  is  equal  to  1 /sampling  frequeney  {f^. 

Notiee  that  Eq.  13  does  not  require  time  derivatives,  so  it  may  handle  noise  better.  These  two 
model  struetures  will  be  eonsidered  and  eompared  in  the  Results  and  Diseussion  seetion. 


The  estimators  eonsidered  here  are  reeursive  least  squares  (RES)  with  exponential  forgetting 
faetor  and  Kalman  filtering  (KE).  Other  methods  whieh  require  a  presumed  ratio  of  time 
constant  will  not  be  used.  The  RES  estimator  assumes  white  noise  (v)  in  the  measurement  in 


The  term  in  the  braeket  is  the  residual  of  the  estimation.  This  is  to  be  weighted  by  a  forgetting 
faetor  X.  Traeking  is  attained  by  plaeing  more  importanee  on  the  more  reeent  data  than  the 
past  data.  A  data  point  is  exponentially  diseounted  as  time  inereases  with  the  forgetting  factor 
(2)  defining  the  time  eonstant  of  the  exponential  deeay.  As  far  as  the  estimator  eoneerns,  there 
are  roughly  1/(1-  2)  of  data  points  whieh  are  effective  (in  memory).  A  lower  value  gives  a 


more  rapid  response  in  tracking  the  parameter  vector.  The  value  commonly  used  in  industry  is 
between  0.97  and  0.995.  In  the  limit  of =  1,  all  data  points  are  remembered,  so  the  estimator 
gives  a  constant  parameter  estimation. 


It  can  be  showed  that  the  following  recursive  formulae  estimate  the  parameter  vector  and 
minimise  the  cost  function: 


where: 

e  is  the  output  (apriori)  prediction  error 
G  is  time  varying  gain  matrix 
P  is  the  covariance  of  parameter  estimation  error 
I  is  the  identity  matrix 

A  is  the  forgetting  factor  with  a  limit  between  zero  and  one 

At  the  start  of  the  estimation,  the  parameter  vector  is  initialized  to  a  zero  vector  and  the 
covariance  P(0)  is  assumed  to  be  10"^I,  which  is  defined  to  be  sufficiently  large,  so  that  fast 
convergence  of  the  solution  can  be  achieved. 

In  a  Kalman  filter,  both  measurement  and  the  parameter  vector  have  white  noise  with 
variances  of  cr^(k)  and  W(k)  respectively,  which  may  change  with  time  step.  Note  that  there 
are  two  parameters  in  the  parameter  vector,  so  W(k)  is  a  2x2  matrix.  The  recursive  formulae 
for  the  Kalman  filter  are  similar  to  RLS,  except  the  gain  and  covariance  matrices  are 
evaluated  according  to  the  noise  variance  as  shown  below. 


where: 


R  is  the  covariance  of  the  parameter  vector. 

al  is  the  estimated  prediction  variance.  It  is  the  sum  of  measurement  noise  variance 
(cr^)  and  prediction  variance  due  to  uncertainty  in  the  parameters  ( cr^  )  as  defined  in 


Eq.  18. 


The  noise  parameters  are  not  assigned  arbitrarily.  They  are  estimated  using  Jazwinski’s 
method  [18].  The  state  noise  covariance  matrix  is  assumed  to  be  an  isotropic  matrix,  W  =  ql. 
The  prediction  variance  (cr^o)  calculated,  assuming  no  state  noise,  i.e.  ^  =  0  by: 


2  2 

If  cr^o  is  smaller  than  the  predietion  error  (e  ),  the  state  noise  is  inferred  to  be  non-zero.  This 
logie  ean  be  implemented  using  a  ramp  funetion  in  Eg,  24, _ 


h(x)  =  • 

x>0 

[o...  otherwise 

20) 
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fire,  q  can  be  updated  according  to: 

-s; 
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Eq.  1 8  shows  that  the  estimated  predietion  varianee  is  made  up  of  the  measurement  noise  and 
the  eomponent  due  to  the  uneertainty  in  the  latent  variables.  Thus  the  measurement  noise  ean 
be  estimated  by  subtracting  the  seeond  eomponent  from  the  prediction  error  as _ 


a\k  +  \)  =  h(s\k  +  \)- x{k  +  \)R{k)x^  {k  + 1)) 

22) 

It  can 

be  proved  that  the  observation  noise  estimated  by  Eq.  22  will  maximize  the  likelihood 

of  the  observation  in  the  next  sample. 

The  Kalman  fdter  is  initialised  with  ^(0)  =  1,  9{0)  =  [0.5,  0.5],  cr^(O)  based  on  40  dB  signal-to- 
noise  ratio.  To  estimate  average  time  constants,  system  noise  is  zero  always  and  is  fixed  to 
have  a  40  dB  signal-to-noise  ratio. 

There  are  two  model  structures  and  two  estimators  forming  a  combination  of  4  techniques. 
These  techniques  will  be  named  according  to  the  model  structure  and  the  estimator  used. 


A  CT  prefix  will  be  used  for  any  technique  that  uses  the  differential  equations  directly  (Eq. 
12).  Techniques  that  use  the  difference  equation  will  have  a  ZOH  prefix.  Kalman  filter  and 
recursive  least  squares  estimators  will  be  abbreviated  as  KE  and  ES  respectively,  and  they 
form  the  suffix  of  a  technique.  Eor  example,  a  CT-KE  technique  is  a  technique  that  uses  the 
differential  equations  and  estimates  the  time  constants  using  Kalman  filtering. 


1,4,4  Temperature  reconstruction 

The  gas  temperature  is  reconstructed  by  applying  Eq.  1  with  the  estimated  time  constants  ( fj 
and  ).  To  make  a  distinction  between  the  actual  gas  temperature  (Eg)  and  the  reconstructed 
temperature,  a  new  symbol  R  is  assigned  to  the  reconstructed  temperature  in  Eq.  23. 


Ei(0  =  7]+fi 


(t)  =  Ej  + 


dm 

dt 

dm 

dt 
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The  assumption  for  two  thermocouple  techniques  is  that  the  two  reconstructed  temperatures 
should  be  equal,  but  not  at  every  time  step.  Noise  and  other  extraneous  factors  may  prevent 
them  being  so.  Most  estimators  assume  that  the  difference  is  a  white  noise.  By  studying  this 
difference  (known  as  residual  of  the  estimation),  this  assumption  can  be  validated.  Any 
departure  from  this  assumption  may  indicate  modelling  error,  biased  time  constant  estimates, 
and  unsteadiness  in  the  thermocouples.  These  issues  will  be  addressed  in  the  model  validation 
section. 


1.5  Results  and  discussion 


1.5.1  Preliminary  evaluation  of  the  data 

Kar  et  al  [11]  pointed  out  that  time  eonstant  estimation  ean  be  eompromised  by  any  DC 
offsets  in  the  raw  thermoeouple  voltages.  This  is  why  the  amplifier  zero  voltages  were 
monitored  and  the  thermoeouple  was  ealibrated  with  the  zero  voltage  taken  off  in  advanee.  An 
aeid  test  is  to  eompare  the  eyelie  means  of  the  two  thermoeouple  signals  to  see  if  there  are  any 
signifieant  DC  offsets,  or  any  other  unaeeounted  effeets  that  violate  the  two  thermoeouple 
assumptions.  For  example,  a  very  steep  temperature  gradient  near  the  thermoeouple  probes 
will  eause  the  two  means  to  differ  from  eaeh  other.  If  so,  the  two-thermoeouple  teehniques 
eannot  be  applied,  and  thermoeouple  eompensation  is  not  possible. 

The  results  were  very  positive.  The  eyelie  mean  of  the  two  thermoeouple  temperatures  are 
differ  by  0.4%  for  helium  (typieal)  and  0.6%  for  nitrogen.  In  one  ease,  the  differenee  is  as  low 
as  0.01%.  The  eyelie  mean  of  TC25  ean  be  bigger  or  small  than  the  one  of  TC50  in  the 
nitrogen  data.  But  the  helium  results  have  shown  that  the  eyelie  mean  of  TC50  is  always 
bigger.  Is  this  systematie  differenee  genuine,  or  this  is  just  a  systematie  error  in  the 
experiment?  In  some  tests,  espeeially  large  stroke  at  high  pressure,  the  zero  readings  were  not 
stable,  so  some  judgment  was  made  on  what  were  the  best  values.  Finally  the  eomparison  is 
not  entirely  foolproof.  If  the  gas  temperature  is  highly  asymmetrie  -  asymmetry  here  refers  to 
dissimilar  waveforms  -  above  and  below  the  mean  level,  it  is  possible  for  the  attenuated 
thermoeouple  signals  with  different  response  times  to  have  different  means.  This  is  beeause 
the  attenuation  is  proportional  to  the  rate  of  ehange  of  the  signal,  not  the  absolute  signal  level. 

1.5.2  Estimated  time  constants 

The  time  eonstants  were  estimated  using  the  four  teehniques  deseribed,  with  and  without  low- 
pass  filtering.  Two  helium  results  at  5  bar  filling  pressures  are  shown  in  Figures  5  and  6.  The 
results  at  8.6  Hz  all  show  smaller  time  eonstants  when  low-pass  filtering  is  not  used.  Many 
previous  studies  [7,  9,  19]  have  found  that  the  effeet  of  measurement  noise  is  to  introduee  a 
bias,  whieh  tends  to  underestimate  the  time  eonstants.  Low-pass  filtering  reduees  the  noise 
level,  so  the  bias  is  redueed  aeeordingly. 

Different  teehniques  handle  the  measurement  noise  differently.  Some  will  give  a  lower  bias 
for  the  same  noise  level.  Sinee  low-pass  filtering  ean  reduee  bias,  the  inerease  in  the  time 
eonstant  after  filtering  ean  reveal  how  well  a  teehnique  handle  noise.  The  results  in  Figure  5 
show  that  the  ZOH-LS  and  ZOH-KF  teehniques  are  more  tolerant  to  noise  than  their 
continuous  time  counterparts.  It  seems  the  model  structure  is  a  more  important  factor  than  the 
estimator  in  noise  rejection. 

The  time  constant  estimates  at  36  Hz  cyclic  frequency  (Figure  6)  are  more  consistent.  The 
ZOH-LS  and  ZOH-KF  both  return  the  same  estimates  of  5.6  ms  (within  0.1  ms)  with  and 
without  filtering.  The  low-pass  fdtering  does  reduce  the  bias  in  Tagawa’s  technique  and  CT- 
LS  technique,  bringing  the  estimates  to  5.6  ms  as  well.  However,  the  increases  in  the  time 
constants  is  very  small  (4%),  which  means  the  bias  is  small  in  the  unfiltered  case.  In 
summary,  the  time  constants  at  36  Hz  are  unaffected  by  bias.  All  techniques  give  similar 
results  and  low-pass  filtering  makes  little  difference.  It  is  of  practical  interests  to  find  out  why 
the  estimation  at  36-Hz  is  of  high  quality,  as  this  may  shed  light  on  how  to  optimise  the 
estimation  process  in  other  less  ideal  situations. 
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Figure  6  Comparison  of  different  two-thermocouple  techniques.  Helium  gas,  5  har  fill  pressure,  10  mm 
stroke  at  36  Hz. 


The  graphs  in  Figure  7  and  Figure  8  can  answer  why  the  estimates  at  36  Hz  have  small  biases. 
The  noise  is  imperceptible  in  the  ensemble  averaged  temperatures.  However,  the  two- 


thermocouple  techniques  either  make  use  of  the  time  derivatives  or  the  temperature  difference 
between  two  time  steps.  Hence,  noise  in  the  time  derivatives  should  be  a  better  measure  of  the 
noise  level.  It  is  clear  that  the  time  derivatives  at  36  Hz  are  much  smoother,  whilst  at  8.6  Hz 
the  time  derivatives  have  lots  of  high  frequency  fluctuations,  which  are  not  physical. 
Obviously,  the  estimators  will  be  more  affected  by  noise  at  8.6  Hz,  i.e.  subject  to  a  lower 
signal-to-noise  ratio,  so  yielding  bigger  biases  as  indicated  by  the  difference  between  the 
estimates  with  and  without  filtering. 
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Figure  7  The  ensemble  averaged  TC25  and  TC50 
temperatures  and  their  time  derivatives.  Helium 
gas,  5  bar  fill  pressure,  10  mm  stroke  at  8,6  Hz, 


Figure  8  The  ensemble  averaged  TC25  and  TC50 
temperatures  and  their  time  derivatives.  Helium 
gas,  5  bar  fill  pressure,  10  mm  stroke  at  36  Hz, 


There  are  several  differences  that  can  explain  why  the  signal-to-noise  ratio  is  higher  at  36  Hz: 

1 .  The  sampling  frequency  and  the  total  sampling  time  are  fixed,  so  if  the  cycle  frequency 
goes  up,  the  number  of  samples  per  cycle  goes  down  with  a  proportional  increase  in  the 
number  of  cycles  sampled.  At  36  Hz,  there  are  360  cycles;  each  one  has  138  samples 
versus  582  samples  for  each  85  cycle. 

2.  Consider  the  time  constants  for  TC25,  the  time  constant  is  about  12  ms  at  8.5  Hz,  which 
means  the  comer  frequency  is  13  Hz.  At  a  cyclic  frequency  of  36  Hz,  the  time  constant 
decreases  to  5.7  ms,  which  means  the  comer  frequency  is  now  28  Hz.  The  cyclic 
frequency  is  relatively  closer  to  the  comer  frequency,  so  less  attenuation  and  stronger 
signals  are  expected. 

3.  Recall  that  the  signals  are  ensemble  averages.  There  are  about  4  times  more  cycles  at  36 
Hz  than  there  at  8.6  Hz,  so  the  signal-to-noise  ratio  should  be  double. 

Amongst  the  above  three  points,  it  has  been  found  that  the  first  point  is  the  most  significant 
factor.  When  any  finite  difference  schemes  calculate  time  derivatives  there  needs  to  be  a  time 
step  {h).  For  example  at  the  time  step  k,  a  forward  difference  scheme  calculates: 


&T{k)  ^  T{k  +  \)-T{k) 

At  h 

24) 

Eq.  24 
noise 

states  that  the  time  derivative  is  inversely  proportional  to  /?.  In  reality,  there  is  also 
[e)  present  in  T.  It  will  be  at  least  as  big  as  in  the  numerator  in  Eq.  24.  Thus, 

AT{k)  ^T{k  +  \)-T{k)  e 

At  h  h 

25) 

Therefore,  any  unwanted  noise  will  also  be  scaled  by  h  in  the  time  derivative.  If  h  is  increased 
by  lower  the  sampling  rate,  the  noise  level  should  be  reduced  accordingly.  The  sampling  rate 
may  be  reduced  by  decimation.  Decimation  by  a  factor  of  r  does  not  simply  mean  selecting 
every  r*^  sample.  This  is  because  doing  so  will  introduce  aliasing  (signals  below  the  Nyquist 
frequency  in  the  original  sampling  rate  is  now  above  the  new  Nyquist  frequency).  To  avoid 
aliasing,  an  anti-aliasing  filter  has  to  be  applied  as  well. 

The  sampling  rate  of  8.5  Hz  data  is  decimated  by  a  factor  of  4  by  firstly  applying  an  8^’’  order 
Chebyshev  filter  with  a  cut-off  frequency  of  0.4  of  the  reduced  sampling  frequency  (0.4  x 
5000/4  =  500Hz).  Only  every  4*^  data  point  in  the  filtered  signals  is  retained.  The  low  pass 
filter  does  not  alter  the  phase  of  the  signal,  and  the  magnitude  is  at  least  99.43%  accurate.  The 
time  derivatives  are  recalculated  and  plotted  in  Figure  9.  Indeed  the  time  derivatives  have 
much  lower  noise,  and  correspondingly  the  biases  in  the  time  constants  are  reduced  or 
eliminated.  As  shown  in  Figure  10,  a  marked  reduction  in  bias  can  be  achieved  just  by  halving 
the  sampling  frequency.  This  happens  for  the  ZOH  model  structure,  even  though  the  time 
derivatives  are  not  used  directly.  As  the  decimation  factor  increases,  the  estimate  eventually 
reaches  a  stable  value,  which  is  identical  to  the  value  found  by  filtering.  In  addition,  the  ZOH 
model  structure  being  more  noise  tolerant,  settles  to  the  stable  value  quicker  than  the 
continuous  time  model  structure  (Tagawa  and  Ohta’s  method).  To  be  consistent,  both  model 
structures  are  processed  using  a  least  squares  estimator,  but  the  influence  of  the  model 
structure  does  not  change  when  a  Kalman  filter  is  used. 
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Figure  9  The  time  derivatives  of  TC25  and  TC50  temperatures,  decimated  by  4  times.  The  sampling 
frequency  is  1.25  kHz.  Helium  gas,  5  bar  fill  pressure,  10  mm  stroke  at  8.6  Hz. 


Decimation  was  applied  to  the  18-Hz  and  36-Hz  data.  The  minimum  decimator  factors 
required  to  achieve  a  stable  time  constants  were  recorded.  It  was  found  that  the  decimation 
factors  were  inversely  proportional  to  the  cyclic  frequency.  For  example,  starting  with  a 
decimator  factor  of  1  at  36  Flz,  it  becomes  2  at  18  Flz,  and  3  at  8.5  Hz.  In  other  words,  the 
sampling  frequency  increases  in  proportion  with  the  cyclic  frequency.  These  results  suggest 
that  the  estimators  can  give  unbiased  estimations  when  the  data  points  per  cycle  are  equal  to 
or  smaller  than  certain  values.  Whilst  it  has  been  shown  that  decimation  can  increase  the 
effective  signal-to-noise  ratio,  it  is  desirable  to  keep  decimation  to  a  minimum  (sampling 
frequency  as  high  as  possible),  so  that  no  genuine  high  frequency  information  is  discarded. 
Also,  the  numerical  error  in  calculating  the  time  derivatives  is  in  the  order  of  h^.  Using  too 
high  a  decimation  will  eventually  lead  to  inaccurate  evaluation  of  the  time  derivatives,  and 
then  the  time  constant  estimates  may  diverge  from  the  stable  values  as  shown  in  Figure  11. 


Figure  10  Time  constant  estimates  of  TC25  as  a  function  of  decimation  factor.  Helium  gas,  5  bar  fill 
pressure,  10  mm  stroke  at  8.6  Hz. 
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Figure  11  Time  constant  estimates  of  TC25  as  a  function  of  decimation  factor.  Helium 
gas,  5  bar  fill  pressure,  10  mm  stroke  at  36  Hz, 

For  each  test  condition,  the  time  constants  were  estimated  with  from  no  decimation  to  10 
times  decimation.  The  minimum  decimator  factor  was  noted  when  the  estimates  had  settled  to 
stable  values.  The  helium  results  are  summarised  in  Table  1.  Only  the  ZOH-LS  results  are 
shown,  since  previous  results  have  established  that  the  ZOFl  structure  handles  noise  better. 
The  estimates  of  Kalman  filtering  (ZOH-KF)  are  within  the  uncertainty  of  the  ZOH  ones.  All 
time  constants  settle  within  0.1  ms  when  the  decimation  increases  further  than  the  stated 
values.  The  results  show  that  the  time  constants  decrease  when  the  filling  pressure  increases 
and  the  cyclic  frequency  increases.  The  cyclic  frequency  seems  to  have  a  bigger  effect. 

Table  1  Time  constant  estimates  nsing  ZOH-LS  techniqne  in  a  helinm  environment.  The  nncertainty 


is  +/-0.1  ms. 


Fill  pressure 

Cyclic 

Time  constant 

Time  constant 

Decimation 

(bar) 

frequency  (Hz) 

of  TC25  (ms) 

of  TC50  (ms) 

factor  used 

1 

8.6 

13.2 

28.0 

5 

5 

8.6 

12.7 

22.8 

4 

18 

8.6 

7.4 

15.4 

2 

1 

18 

8.6 

22.1 

2 

5 

18 

7.3 

17.1 

2 

18 

18 

6.2 

15.0 

3 

1 

36 

8.4 

23.6 

3 

5 

36 

5.6 

14.8 

2 

18 

36 

5.7 

14.5 

2 

Table  2  gives  the  nitrogen  results.  The  magnitude  of  the  time  constants  is  3  to  4  times  bigger 
than  that  of  helium.  So,  the  raw  temperature  is  a  lot  more  attenuated.  Probably  due  to  the 
increases  in  time  constant  magnitude,  there  are  more  variations  in  the  time  constant  estimates 
as  the  decimation  factor  increases.  The  uncertainty  in  Table  2  refers  to  the  possible  change  in 
the  estimates  as  the  decimation  factor  varies.  The  response  of  time  constants  to  fill  pressure 
and  frequency  is  not  as  clear  as  with  helium.  In  particular,  the  time  constants  at  5  bar  and  36 
Hz  are  exceptionally  high.  The  uncertainty  is  very  high  too. 


Table  2  Time  constant  estimates  nsing  ZOH/LS 

techniqne  in  a  nitro 

gen  environment. 

Fill 

pressure 

(bar) 

Cyclic 

frequency 

(Hz) 

Time  constant 
of  TC25  (ms) 

Time 

constant  of 
TC50  (ms) 

Uncertainty 

(ms) 

Decimation 
factor  used 

1 

8.6 

46.0 

119.3 

+0.0, -0.4 

3 

5 

8.6 

37.0 

95.1 

+0.1, -0.1 

3 

18 

8.6 

40.8 

95.3 

+0.0, -0.1 

3 

1 

18 

39.4 

110 

+0.2,  -0.4 

4 

5 

18 

37.2 

101.8 

+0.1,  -0.8 

3 

18 

18 

28.3 

73.3 

-0.0,  -0.8 

3 

1 

36 

40.2 

117 

+0.0,  -2.0 

2 

5 

36 

62 

170 

+ 

o 

■kl 

2 

18 

36 

24.9 

59.8 

+0.0,  -1.7 

1 

1,5,5  Reconstructed  temperatures 

Ultimately  the  time  constant  is  used  for  finding  the  true  gas  temperature.  The  reconstructed 
temperatures  are  symbolised  as  R25  for  TC25  and  R50  for  TC50.  This  is  to  avoid  confusion 
with  the  actual  gas  temperature  at  the  regenerator  exit.  The  estimated  values  are  listed  in 
Table  1  and  Table  2.  The  time  derivatives  are  calculated  using  central  difference  scheme. 

The  reconstructed  temperatures  at  three  frequencies  are  plotted  in  Figures  13-15.  Results  at 
other  conditions  are  given  at  the  end  of  this  report.  The  results  for  R25  and  R50  are  very 
similar  (Figure  12),  so  only  R25  are  shown.  Although  the  true  gas  temperature  is  not 
available,  the  accuracy  of  the  reconstructed  temperatures  can  be  evaluated  by  the  common 
features  in  these  three  plots.  Firstly,  the  two  techniques  ZOH-LS  and  ZOH-KF  give  almost 
identical  temperatures  and  this  is  an  assurance  of  the  quality  of  the  reconstruction.  Secondly, 
the  peak  reconstructed  temperature  is  in  phase  with  the  cold  wire  temperature,  which  indicates 
the  phase  of  the  thermocouple  signal  has  been  corrected  properly.  Phase  correctness  is  only 
true  if  the  regenerator  has  little  effect,  i.e.  the  temperature  changes  are  driven  by  compression. 
The  reconstruction  and  cold  wire  temperatures  are  in  phase  in  the  first  quarter  of  the  cycle. 
This  is  when  the  piston  is  moving  up,  and  gas  mixing  in  compression  space  is  not  expected, 
so  the  regenerator  cannot  affect  the  cold  wire  temperature  much. 


Figure  12  Reconstructed  temperatures  using  ZOH-LS  technique.  Helium  gas,  5  bar  fill  pressure,  10  mm 
stroke  at  8.6  Hz. 
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Figure  13  Comparison  of  the  raw  and  compensated  thermocouple  and  cold  wire  temperatures.  The 

sampling  frequency  was  decimated  to  1.25  kHz.  Helium  gas,  5  bar  fill  pressure,  10  mm  stroke 
at  8.6  Hz. 
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4  Comparison  of  the  raw  and  compensated  thermocouple  and  cold  wire  temperatures.  The 
sampling  frequency  was  decimated  to  2.5  kHz.  Helium  gas,  5  har  fill  pressure,  10  mm  stroke  at 
18  Hz. 


Figure  15  Comparison  of  the  raw  and  compensated  thermocouple  and  cold  wire  temperatures.  The 
sampling  frequency  was  2.5  kHz.  Helium  gas,  5  har  fill  pressure,  10  mm  stroke  at  36  Hz. 


1.6  Model  validation 

The  thermocouple  measurements  are  compensated  using  the  estimated  time  constant.  In  other 
words,  these  compensated  signals  are  no  longer  real  measurements.  They  are  only  valuable  if 
they  can  be  validated  by  some  means.  The  obvious  method  is  to  compare  them  against  the 
true  temperature,  but  this  is  clearly  impossible.  In  practice,  the  compensated  signals  can  be 
indirectly  validated  by  model  validation.  Model  validation  checks  that  the  thermocouple 
models  are  sufficient  to  describe  the  response  of  the  thermocouple  and  other  assumptions 
made  in  the  estimators.  The  important  assumptions  made  in  two-thermocouple  techniques  are 
summarised  below: 

•  Convective  heat  transfer  is  dominant,  such  that  a  first  order  response  equation  is 
sufficiently  accurate. 

•  The  two  thermocouples  measure  the  same  temperature  field,  i.e.  the  separation  distance 
between  them  is  small  enough. 

•  Once  the  time  constant  estimate  converges  as  the  decimation  factor  increases,  the  estimate 
is  optimal  (least  biased  or  unbiased). 

•  Measurement  noise  and  the  differences  in  the  two  reconstructed  temperatures  are  due  to 
white  noise. 

•  The  average  time  constants  are  sufficiently  good  that  time  varying  techniques  are  not 
necessary. 


The  existence  of  noise  means  that  and  R2  cannot  be  equalled  at  all  time.  If  all  the 

assumptions  hold,  the  residual  in  the  estimation  will  be  non  zero,  and  randomly  distributed 
with  zero  mean,  i.e.  a  white  noise.  Residual  {aposterior  prediction  error)  is  defined  as 


ri{k)  -  y(k)  -  (k)d(k) 


26) 


Model  validation  is  essentially  a  white  noise  test.  It  finds  the  autocorrelation  function  of  the 
residual,  .  The  autocorrelation  function  compares  if  the  present  sample  (zero  delay)  to 

sample  at  a  later  time  step  (non-zero  delay).  If  the  residual  is  white  noise,  it  should  be 
uncorrelated  to  anything  else  except  itself  So,  the  autocorrelation  function  whose  value  is 
normalised  to  the  value  at  zero  delay  should  be  one  at  zero  delay  and  zero  everywhere  else.  In 
practice  however,  because  the  data  sample  is  finite,  this  limits  the  noise  bandwidth.  The 
consequence  is  that  the  samples  may  be  slightly  correlated.  A  95%  confidence  interval  has 
been  constructed  to  allow  for  the  effect  of  finite  data  length.  If  the  autocorrelation  exceeds 
this  interval,  this  means  it  is  95%  certain  that  the  residual  is  correlated  (not  white). 

It  has  been  deduced  that  when  no  decimation  or  low-pass  filtering  is  applied,  the  time 
constant  estimates  tend  to  be  biased  due  to  the  influence  of  noise.  Figure  16  and  Figure  17 
compare  the  autocorrelation  function  of  residual  without  and  with  decimation.  The  decimation 
selected  was  the  minimum  decimation  required  to  achieve  a  stable  time  constants.  The  results 
clearly  show  that  when  the  time  constant  was  too  low,  there  are  dynamics  in  the  residual 
(unaccounted  system  information)  which  are  correlated.  The  function  is  always  above  the 
upper  limit  of  uncertainty.  After  appropriate  decimation  is  applied,  the  autocorrelation 
fluctuates  around  zero  (except  delay  =  0),  and  the  fluctuations  are  within  the  confidence 
interval.  The  residual  can  be  accepted  as  white  noise.  This  test  also  indicates  that  the 
thermocouple  model  is  sufficiently  accurate  and  that  any  effect  of  thermocouple  separation  is 
negligible. 
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Figure  16  Normalised  Autocorrelation  function  of  the  residual  for  the  CT-LS  technique.  Nitrogen  gas,  1 
har  fill  pressure,  10  mm  stroke  at  8.6  Hz.  No  decimation  was  applied.  The  dotted  lines  show 
the  95%  confidence  interval  within  which  the  autocorrelation  should  lie  if  the  residual  is  white 
noise. 
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Figure  17  Normalised  autocorrelation  function  of  the  residual  for  the  CT-LS  technique.  Nitrogen  gas,  1 
har  fill  pressure,  10  mm  stroke  at  8.6  Hz.  8  times  decimation  was  applied.  The  dotted  lines 
show  the  95%  confidence  interval  within  which  the  autocorrelation  should  lie  if  the  residual  is 
white  noise. 


The  results  of  the  model  validation  are  generally  satisfaetory  at  higher  frequeney  and  filling 
pressure,  see  Figure  18  for  example.  But,  there  are  some  marginal  eases  where  the 
autoeorrelation  funetion  exeeeds  the  eonfidenee  intervals  by  a  small  amount.  The  ease  for 
nitrogen  gas  at  18  bar  filling  pressure  at  8.6  Hz  elearly  illustrates  this  point  (Figure  19).  The 
marginal  ease  should  still  be  eonsidered  aceeptable  beeause  there  is  a  5%  probability  that  a 
genuine  white  noise  fails  the  model  validation.  If  the  marginal  tests  are  aeeepted,  8  out  of  9 
tests  pass  the  model  validation  when  nitrogen  is  used,  but  only  3  out  of  9  tests  pass  for 
helium.  Table  3  summarises  the  results  for  all  tests.  Amongst  all  the  failed  eases,  there  is  a 
eommon  feature  -  the  autoeorrelation  funetion  only  fails  at  delay  =  3.  For  all  other  delays  the 
funetion  is  within  the  limits.  An  example  of  this  failure  is  shown  in  Figure  20.  The  failure  at 
delay  =  3  means  that  the  model  does  not  aeeount  for  the  dynamies  with  a  delay  of  3.  The 
results  eannot  pinpoint  whieh  assumption  has  been  violated.  For  these  eases,  variable  time 
eonstants  were  estimated  instead  of  average  values,  the  autoeorrelation  does  not  ehange  much. 
This  shows  that  the  high  correlation  at  delay  =3  is  not  caused  by  the  average  time  constant. 
Grys  et  al.  [15]  performed  similar  model  validation,  in  which  a  least  mean  squares  estimator 
was  used.  Interestingly  the  normalised  autocorrelation  function  also  has  a  distinct  spike,  but  at 
delay  =  2  (Figure  21).  Nevertheless  Grys  et  al.  [15]  still  stated  that  “the  autocorrelation  lies 
inside  the  confidence  interval”.  It  is  not  apparent  at  this  stage  whether  this  feature  in  the 
autocorrelation  is  enough  to  jeopardise  the  two-thermocouple  technique.  Since  the  actual  gas 
temperature  is  unknown,  it  is  difficult  to  establish  which  other  assumptions  are  to  blame,  and 
decide  appropriate  action.  Therefore,  it  is  recommended  to  proceed  with  caution.  Without 
rigorous  analysis,  it  seems  likely  that  the  conduction  effect  may  play  a  role  which  the 
thermocouple  model  needs  to  take  account  of 

Table  3  Results  of  the  model  validation  tests  using  average  time  constant  estimation,  ZOH-LS 


techniques  results 


Filling 

pressure 

(bar) 

Cyclic 

frequency 

(Hz) 

Helium  results 

Nitrogen  results 

1 

8.6 

Delay  =  3  exceeds  the  limit 

Acceptable 

5 

8.6 

Delay  =  3  exceeds  the  limit 

Marginal 

18 

8.6 

Delay  =  2  and  3  exceed  the  limit 

Marginal 

1 

18 

Failed.  The  autocorrelation 
function  appears  to  be  periodic 

Acceptable 

5 

18 

Delay  =  3  exceeds  the  limit 

Delay  =  3  exceeds  the  limit 

18 

18 

Marginal 

Marginal 

1 

36 

Marginal 

Acceptable 

5 

36 

Acceptable 

Acceptable 

18 

36 

Delay  =  3  exceeds  the  limit 

Acceptable 

1 


Delay 

Figure  18  Normalised  autocorrelation  function  of  the  residual  using  ZOH-LS  techniques.  No  decimation 
has  been  applied.  Nitrogen  gas,  18  har  fill  pressure,  10  mm  stroke  at  36  Hz.  The  dotted  lines 
show  the  95%  confidence  interval  within  which  the  autocorrelation  should  lie  if  the  residual  is 
white  noise. 


Delay 


Figure  19  Normalised  autocorrelation  function  of  the  residual  using  ZOH-LS  techniques.  3  times 

decimation  has  been  applied.  Nitrogen  gas,  18  bar  fill  pressure,  8  mm  stroke  at  8.6  Hz.  The 
dotted  lines  show  the  95%  confidence  interval  within  which  the  autocorrelation  should  lie  if 
the  residual  is  white  noise. 


1 


Delay 

Figure  20  Normalised  autocorrelation  function  of  the  residual  using  ZOH-LS  techniques.  2  times 

decimation  has  been  applied.  Nitrogen  gas,  5  har  fill  pressure,  10  mm  stroke  at  18  Hz.  The 
dotted  hues  show  the  95%  confidence  interval  within  which  the  autocorrelation  should  lie  if 
the  residual  is  white  noise. 
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Figure  21  Normalised  autocorrelation  function  of  the  residual  from  Gry  et  al.  [15].  k  represents  delay 
and  a  is  the  error  prohahility.  1-a  indicates  the  width  of  the  confidence  interval. 


The  helium  result  at  1  bar  pressure  and  1 8  Hz  deserves  further  discussion  because  it  is  the  one 
that  fails  differently.  The  autocorrelation  function  shows  regular  oscillations  in  Figure  22.  The 
confidence  interval  is  consistently  exceeded  over  a  range  of  delays.  This  is  in  a  sharp  contrast 
to  the  case  in  which  the  autocorrelation  is  unusually  high  at  delay  =  3.  Therefore,  the  time 
constants  and  the  reconstructed  temperature  at  this  condition  may  be  unreliable  and  are 
ignored  accordingly. 


Figure  22  Normalised  autocorrelation  function  of  the  residual  using  ZOH-LS  techniques.  2  times 
decimation  has  been  applied.  Helium  gas,  1  har  fill  pressure,  10  mm  stroke  at  18  Hz.  The 
dotted  lines  show  the  95%  confidence  interval  within  which  the  autocorrelation  should  lie  if 
the  residual  is  white  noise. 


1,7  Conclusions 

The  following  conclusions  can  be  drawn  from  the  thermocouple  compensation  work: 

•  Although  very  fine  thermocouples  were  used  in  the  experiment,  thermocouples  of 
different  diameters  still  had  different  responses,  which  suggested  the  time  response 
was  still  insufficient  to  measure  the  gas  temperature  accurately. 

•  The  cyclic  variability  is  found  to  be  low.  This  makes  ensemble  averaging  possible 
without  losing  much  information. 

•  Low-pass  filtering  reduces  bias  in  time  constant  estimation,  so  comparing  the  results 
with  and  without  low-pass  filtering  can  reveal  the  noise  tolerance  of  different 
techniques.  The  results  show  that  model  structure  is  more  important  in  noise  rejection 
than  the  types  of  estimator  used.  The  ZOH  model  structure  is  far  superior  to  the 
continuous  time  (CT)  model  structure. 

•  Low  noise  as  found  in  high  cyclic  frequency  results  has  prompted  the  fact  that  any 
unwanted  noise  will  be  scaled  by  the  time  step  in  the  time  derivation.  So,  an  effective 
way  to  reduce  noise  is  to  reduce  the  sampling  frequency  by  decimation. 

•  It  has  been  found  that  the  time  constant  estimates  converge  as  the  decimation  factor 
increases.  The  stable  value  allows  a  minimum  decimation  factor  to  be  determined, 
above  which  an  unbiased  time  constant  estimate  can  be  obtained. 

•  The  time  constants  are  found  to  decrease  when  the  filling  pressure  increases  and  the 
cyclic  frequency  increases.  The  cyclic  frequency  seems  to  have  a  bigger  effect. 

•  The  magnitude  of  the  time  constants  in  a  nitrogen  environment  is  3  to  4  times  bigger 
than  that  in  helium. 


•  Features  in  the  reeonstrueted  temperatures,  e.g.  the  peak  temperature  in  phase  with  the 
eold  wire  temperature,  suggest  that  the  eompensation  is  physieally  sound. 

•  8  out  of  9  nitrogen  results  pass  the  model  validation  test,  but  only  3  out  of  9  helium 
results  pass.  Most  of  the  failure  ean  be  eharaeterised  by  abnormally  high 
autoeorrelation  at  delay  =  3. 

•  Variable  time  eonstants  were  estimated  for  the  results  whieh  fail  the  model  validation 
test.  However,  the  autoeorrelation  does  not  ehange  mueh  when  eompared  to  the 
average  time  eonstant  eases. 
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1,9  Appendix 


#he#1  p1s10f_08.58.xls,  No  filtering 


#he#1p18s10f_08.58.xls,  No  filtering 


Helium  gas,  1  bar  fill  pressure,  1 0  mm  stroke 


Helium  gas,  18  bar  fill  pressure,  10  mm 


at  8.6  Hz. 


#he#1p1s10fJ8.04.xls,  No  filter 


stroke  at  8.6  Hz. 

#he#1p18s10fJ8.03.xls,  No  filter 


Helium  gas,  1  bar  fill  pressure,  1 0  mm  stroke 


Helium  gas,  18  bar  fill  pressure,  10  mm 


at  18  Hz. 


#he#1p1s10f_36.22.xls,  No  filtering 


stroke  at  18  Hz. 

#he#1p18s10f_36.05.xls,  No  filtering 


Helium  gas,  1  bar  fill  pressure,  1 0  mm  stroke 


Helium  gas,  18  bar  fill  pressure,  10  mm 


at  36  Hz. 


stroke  at  36  Hz. 


#n2#1p1s10f_08.62.xls,  No  filter 


#n2#2p18s8f_08.57.xls,  No  filtering 


Nitrogen  gas,  1  bar  fill  pressure,  10  mm 
stroke  at  8.6  Hz. 


Nitrogen  gas,  18  bar  fill  pressure,  8  mm 
stroke  at  8.6  Hz.  10  mm  was  not  achievable. 


#n2#1p1s10fJ8.12.xls,  No  filter 


#n2#1p18s10f_18.35.xls,  No  filter 


Nitrogen  gas,  1  bar  fill  pressure,  10  mm 
stroke  at  18  Hz. 


Nitrogen  gas,  18  bar  fill  pressure,  10  mm 
stroke  at  18  Hz. 


#n2#1p1s10f_36.13.xls,  No  filter 


#n2#1p18s10f_36.40.xls,  No  filter 


Nitrogen  gas,  1  bar  fill  pressure,  10  mm 
stroke  at  36  Hz. 


Nitrogen  gas,  18  bar  fill  pressure,  10  mm 
stroke  at  36  Hz. 


#n2#1  p5s10f_08.57.xls,  No  filtering 


Nitrogen  gas,  5  bar  fill  pressure,  10  mm  stroke  at  8.6  Hz. 


#n2#1p5s10f_18.05.xls,  No  filter 


Nitrogen  gas,  5  bar  till  pressure,  10  mm  stroke  at  18  Hz. 


#n2#1p5s10f_35.93.xls,  No  filter 


Nitrogen  gas,  5  bar  till  pressure,  10  mm  stroke  at  36  Hz. 


